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The influence of continuous association of one component on the free energy of a non- 


electrolyte solution is derived from a few assumptions. The result is in accord with available 
data for hydrocarbon-alcohol solutions. Systems of two associating components and the general 
interpretation of the free energy of non-electrolyte solutions are briefly discussed. 


1. INTRODUCTION 


HE examination of the thermodynamic 

properties of non-electrolyte solutions re- 
veals a variety of characteristic types. The 
interpretation of these types is desirable both for 
the extension of our basic knowledge and for 
technical applications such as the design of 
efficient distillation columns. 

Systems of an associating and a non-associ- 
ating component represent one of these types, 
easily recognized in a diagram of the function 
log(y/yvo) (y and yo activity coefficients) against 
the mole fraction x. The advantages of the 
examination of this particular function and its 
representation by 


log(y/vo) = B(1 — 2x) +C[—1+6x(1—x)] 
(1) 


have been discussed in the second paper of this 
series.! It has been found that associating systems 
can be approximated by (1) with appreciable 
values of the coefficients B and D, while C is 
small. 

! To be published in Ind. and Eng. Chem. 
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It is the object of this paper to derive this 
peculiar behavior from a suitably chosen model. 
In other words, an attempt is made to find the 
relationship between the free energy and the 
stoichiometric mole fraction x for a solution 
which consists of molecules H of a non-associating 
substance like a hydrocarbon and monomeric and 
associated molecules A, As, A3:--A,-:: of an 
associated substance like an alcohol. 

Obviously this problem can be solved only by 
means of an approximation which is based on 
far-reaching simplifying assumptions. The de- 
cisive point is the proper selection of the assump- 
tions. If they are not far-reaching enough, the 
calculation becomes too complicated and the 
number of empirical constants so large that the 
result has no value at all. On the other hand, the 
assumptions should not slur over any charac- 
teristic feature, and above all, they must be 
consistent. 

The relations derived in this paper are based 
on three assumptions, each of which is more or 
less well supported by previous knowledge. They 
are: 
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(a) The free energy of interaction can be 
approximated by Scatchard’s method.!* 

(b) The non-ideal or interaction parts of partial 
molal quantities (free energy, heat content, 
volume) of members of a homologous series are 
proportional to the molecular size. This is a 
somewhat generalized form of a hypothesis which 
has been introduced first by Kratky and Musil? 
in a study of the osmotic pressure of high poly- 
mers. It has been used since by many authors in 
this field. 

The justification of this assumption lies in the 
fact that the non-ideal functions directly repre- 
sent the interaction between molecules, and that 
any part of a molecule (like an OH-group or a 
CHe2-group) interacts with its environment ap- 
proximately the same way irrespective of being 
part of a larger or a smaller molecule. 

(c) It is assumed that the association constant 
for the reaction 


A+Ajy-1=Ay 


does not depend on the order g of the association. 
This hypothesis has been considered by Las- 
settre,* later used by Kempter and Mecke‘* and 
thoroughly discussed by other authors.® The basic 
idea is similar to that of the second assumption. 

The three assumptions differ in significance. 
Assumptions (a) and (b) are fair approximations 
required for a reasonably simple derivation. 
Assumption (c), however, restricts the discussion 
to a certain type of association. Other types 
undoubtedly exist. For instance, for a solution 
of acetic acid in a hydrocarbon association to 
double molecules rather than “continuous” 
association should be expected. 

Starting from different assumptions, Flory** 
and Tobolsky and Blatz®» obtained a different 
result for a related problem. 

1a G. Scatchard, Chem. Rev. 8, 321 (1931); J. Am. Chem. 
Soc. 56, 995 (1934); Trans. Faraday Soc. 33, 160 (1937). 

20. Kratky and A. Musil, Oesterr. Chem. Ztg. 40, 144 
(1937); Z. Elektrochemie 43, 326 (1937). 

3E. N. Lassettre, J. Am. Chem. Soc. 59, 1383 (1937); 
Chem. Rev. 20, 281 (1937). 

4H. Kempter and R. Mecke, Zeits. f. physik. Chemie 
(B) 46, 229 (1940). ; 

5 J. N. Wilson, Chem. Rev. 25, 400 (1939); R. S. Ras- 
mussen, Thesis, California Institute of Technology, Pasa- 
dena (1941); J. Kreuzer and R. Mecke, Zeits. f. physik. 
Chemie (B) 49, 309 (1941); E. G. Hoffman, ibid. 53, 179 
(1943); J. Kreuzer, ibid. 53, 213 (1943). 

5a P. J. Flory, J. Chem. Phys. 14, 49 (1946). 


5» A. V. Tobolsky and P. J. Blatz, J. Chem. Phys. 13, 
379 (1945). 


2. NOTATION 


The present problem requires a clear distinc- 
tion between stoichiometric and “true’’ mole 
numbers and fractions. The symbol x is used for 
the stoichiometric mole fraction of the associ- 
ating substance (alcohol) so that 1—x represents 
the stoichiometric mole fraction of the non-asso- 
ciating component (hydrocabon). The “true” 
mole numbers, etc., are denoted according to the 
following scheme: 


molecule H A_ Az Ages: 
mole fraction No Ni Nez 
activity coefficient Bo 6: Bz Bs By 


We adopt the rule that the symbol >> indicates 
summation over the values 1, 2, 3--- of the 
subscript, not including zero. 


With 
n= yn; N=DN,; N*=)> gN, (2) 


the relations between the quantities so far 
introduced are 


No Nh 
No= ; M= ; N= =1—Np; 
No+n Notn Notn 
Nea g (3) 
Notn 
Lign, 
x = (4) 


mo+Xen, 


since a polymerisate of g single alcohol molecules 
contributes g molecules to the analytical or 
stoichiometric mole number. The last equation 
can be transformed, by means of (3), in the 
following way 


N* No x N* 
Not+tN* 1-x 1-N 


5 
(S) 


The significance of the quantity N and nu- 
merous related functions has been discussed at 
very great length by several authors. Various 
elaborate systems of nomenclature were sug- 
gested. It is obvious that a ratio like N*/N is 
some measure of the degree of association, and 
that other functions can be constructed, any of 
which characterizes the state of association. 
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3. THE “TRUE” MOLE FRACTION 


Since the free energy will be expressed as a 
function of the ‘‘true’’ mole fraction JN, the rela- 
tionship between this function and the stoichio- 
metric mole fraction x is required. It can be 
derived from the assumptions (b) and (c). 

According to (c), the association equilibrium 
is characterized by 


N By = (6) 


where K has the same value for all values of g. 
Repeated application of this equation to the 
values g—1, g—2---2 and multiplication of the 
obtained relations leads to 


N,By= (7) 


According to assumption (b), the true activity 
coefficients are.connected by 


logB, = g (8) 
so that 
N,=KNy. (9) 


Similar formulas have been previously used 
without justification by assumption (b). 
Equations (2) and (9) furnish 


=> (10) 


Ni=N/(1+KN); 


=> (11) 


so that, with the aid of (2) and (9) 


N*¥=)) g(KN1)*/K 
Equation (10) furnishes 


N(1+KN) =Ni/(1—KNj)?. (13) 


From (5), (12), and (13) we obtain 
(14) 
or, arranged in powers of N, 
(15) 
so that 
N={-—1+[1+4Kx(1—x) }}/[2K(1—x)]. (16) 
The negative sign of the root is excluded by the 


fact that N must be positive. This formula 
satisfies, of course, the obvious conditions 


N=x for K=0; (17) 
N=0 for x=0; (18) 
N=1 for x=1. (19) 


The application of this relation will be discussed 
later. 


4. THE FREE ENERGY 


The free energy will be represented in two 
ways. First, in the conventional way by means of 
the stoichiometric and experimental quantities: 


(F— F°)/2.3RT =(1—<x) log(1—x) 
+x logx+Q, (20) 


where the function 


Q=(1—x) logyo+x logy (21) 


’ 


corresponds to Scatchard’s ‘‘excess free energy.” 
The activity coefficients of the hydrocarbon and 
the alcohol are denoted by 7» and y, respectively. 
The quantities F and Q refer to that amount of 
solution which contains (1—x) moles of hydro- 
carbon and x moles of alcohol. 

The free energy can also be represented in 
terms of the “‘true’’ mole fractions: 


(F— F°)/2.3RT =(No logNo 
+> N, logN,+Q*)/(Not+N*), (22) 


Q* = No logBot+ N, logB,, (23) 


where the “‘true”’ activity coefficients are denoted 
by Bo and B,. The denominator No+N* reduces 
the right hand side of (22) to the amount of 
solution containing 1—x and x moles of the 
components since the numerator refers to No 
moles of hydrocarbon and }-gN,=N* moles of 
monomeric alcohol. 

According to assumption (a) the quantity Q* 
can be represented by Scatchard’s equation 


Q* bog Vo bon Va 
(VoNo+d (24) 


where bo, denotes Scatchard’s coefficient of inter- 
action between a solvent molecule and a polymer 
of order g, divided by 2.3RT. The coefficient },, 
indicates the interaction of two polymers of 
order g and h. Since Q* is a homogeneous function 
of first degree in No and Nj, the activity coef- 
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ficients are given by 
logBo=0Q*/9No; logB,=dQ*/AN,. (25) 


It can be concluded from Tyrer’s measure- 
ments® of the densities of methanol-benzene 
solutions that the apparent molal volume of the 
alcohol varies not more than a few tenths of one 
percent over the whole range. The approximation 


is therefore amply sufficient. By introducing (26) 
into (24) and carrying out the differentiation 
indicated in (25), one recognizes that condition 
(8) requires that the interaction coefficients bo, 
and b,, be independent of the subscripts g and h, 
i.e., that 

bog=b01; (27) 


whatever the values of g and h are. 
Equations (24) to (27) furnish 


logBo Vo(b01— 011) 2"; 
logB, = g (bo1—b11) (1 — 2”) +11], 
where 


ViN*/( VoNot+ ViN*) 
=Vix/LVo(l—x)+Vix] (29) 


represents the volume fraction of the alcohol. 
The change from the ‘‘true’”’ to the experi- 
mental or stoichiometric activity coefficients can 
be carried out by means of (20) to (23). More 
convenient is a direct calculation based on the 
fact that the activities, as experimental quan- 


(28) 


TABLE I. Values of the association function A. 
(— sign for 0<x<3$; + sign for }<x<1.) 


K « 0.05 0.1 0.2 0.3 0.4 
Fo 0.95 0.9 0.8 0.7 0.6 


1 0.0107 0.0163 0.0199 0.0165 0.0090 
2 -0290 0433 .0494 -0396 
3 0486 0702 0777 -0330 
-0962 -1036 0429 
6 


-0870 -1206 0976 -0519 
1049 1475 -1126 -0600 


8 -1386 1835 1854 1397 0738 
10 -1693 -2193 -1625 -0856 
12 1974 2512 2457 
14 2233 -2800 -2001 .1049 
16 .2472 3063 2933 -2160 -1130 
20 -2904 -3528 -3328 -2434 -1268 
25 3376 4024 2722 -1416 
30 3786 4450 4094 .2966 1512 


6D. Tyrer, J. Chem. Soc. 99, 871 (1911). 


tities, do not depend on any assumption regard- 
ing the molecular state. Denoting the activities 
in the reference states by do® and a,°, we have, 
therefore, 


xy=NiB:/a1°. (30) 
From (3), (5), (10), and (14), we derive 
N./(1—x) =No+N*=1—N+N* 

=1+KN?*; (31) 
Ni N NotN* 1+KWN? 


x 1+KN N*  (1+KWN)? 


(32) 


Equations (28) to (32) furnish 


logyo=log(1+KN?) 
+ (33) 
logy 
+ (1 —z)?+bi; |+loga,°. (34) 
The reference states are defined by the usual 
convention that 


logyo=0 for x=0; logy=0 for x=1, (35) 


the corresponding limits N=0 and N=1 being 
given by (18) and (19). Therefore, 


logy 
+ Vi(do1 —by)(1 (36) 


The function log(y/yo) is useful for a com- 
parison with experimental data. With the ab- 
breviation 


A=2x log(i+K)—2log(it+KN) (37) 
it is represented by 


log(y/v0) =A +(1—2x) log(1+K) 
—z)?— J. (38) 


Developing the expression in brackets in terms of 
v=(Vi-—Vo)/(Vit Vo), (39) 
one obtains 
Vi(1 —z)? Vo2? Vi Vo/( Vit Vo) 
X (40) 


The term containing v is a second order term 
which is outside the scope of the present dis- 
cussion. Indeed, the error introduced by (27) can 
be estimated in a somewhat complicated calcu- 
lation and is found to be of the same order of 
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magnitude as the v-term in (40). This term, 
therefore, will be neglected or taken care of by 
means of an empirical coefficient so that with 


Vo) (41) 
the result 
log(y/vo) =A+B(1—2x) (42) 


is obtained. 

The linear B-term represents the usual first 
approximation of the interaction between the 
two components. The coefficient B, however, 
contains log(i+<X) and is to be expected, there- 
fore, to be unusually large for associating 
systems. The “‘association function’’ A depends 
on the association constant K and the mole 
fraction x in a complicated way and determines 
the peculiar behavior of systems with an asso- 
ciating component. 


5. THE ASSOCIATION FUNCTION 


The association function (37) is odd with 
respect to the variable (1—2x). In other words, 
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6 
5 
j 
° 0.1 0.2 x 0.3 0.4 0.5 
1.0 0.9 0.8 0.7 0.6 
x<0. fo) 
? 2 


if A’ represents the value of A for x’, and A” 
the value for x”’=1—x’, then A’+A”=0. To 
prove this, we calculate first 
1+KN=(1—2x+7r)/2(1—x); 
r=+[1+4Kx(1—x)]! (43) 


from (16). Then we see that, since r’=r’’, 
(1+KN’)(1+KN”) =(1—2x-+,7) 
(44) 


after an appropriate rearrangement. Equation 
(37) then furnishes 


A’ +A" =2(x+1—x) log(1+K) 
—~2log(1+K)=0. (45) 


For small values of K, the square root 7 can 
be developed into a power series. The develop- 
ment up to the second order furnishes 


A=K?x(1—x)(2x—1)/2.303, (46) 


valid for small values of K and any value of x. 
Likewise one obtains for small values of x and 
any value of K 


A = 2x[log(1+K) —K/2.303]. (47) 


If K<1 and therefore (46) holds, Eq. (42) is 


Fic. 1. Toluene-methanol = n-heptane-methanol 
mm). 


represented by (1) with 
C=0; -—8D=K?/2.303, (48) 


i.e., the association function furnishes the 
D-term. For greater values of K appreciable 
deviations arise though the general course is 
still similar. It has been realized already by 
Scatchard that association leads to a term of 
this type. 

Table I contains values of the function A. 


6. EXPERIMENTAL DATA 


The results of Benedict’ et al. on toluene- 
methanol and mn-heptane-meihanol at atmos- 
pheric pressure are shown in Fig. 1. In this fig- 
ure, the absolute magnitudes of the function (42) 
are plotted against the mole fraction x of the asso- 
ciating component in such a way, that the values 
x and (1—.) of the abscissa coincide. In a diagram 
of this kind, the two branches of [log(y/yo) ] are 
represented by a single curve if, and only if, log 


7™M. Benedict, C. A. Johnson, E. Solomon, and L. C. 
Rubin, Trans. Am. Inst. Chem. Eng. 41, 371 (1945). 
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(y/vo) is an odd function of (1—2vx), i.e., if the 
values of log(y/7o) for x and (1—x) have opposite 
sign and equal magnitude. 

The data of Benedict et al. were checked by 
means of the method described in the first paper 
of this series and found very well consistent. The 
scattering, as shown in Fig. 1, is small. The dia- 
gram shows immediately that the two branches 
of [log(y/vo)] for x«<0.5 and x>0.5 can be 
represented by a single curve which deviates 
considerably from a straight line. The curve for 
toluene-methanol in Fig. 1 corresponds to Eq. 
(42) with B=0.95 and K=3.8 in the asso- 
ciation function (37). This curve represents the 
data well within the small experimental errors. 
The data for n-heptane-methanoi are fairly well 
represented by the curve in Fig. 1 corresponding 
to B=1.25 and K=6.6 

The results of Ferguson® on m-hexane-methanol 
at 45°C can be represented by B=1.71 and 
K=30. Data for n-heptane-ethanol® at 30°C 
scatter considerably more. Even so, the expected 
influence of the association of ethanol is unmis- 
takable. 

In these examples the coefficient C (Eq. (1)) is 
practically zero. Frequently, however, an addi- 
tional, usually small, C-term is required for the re- 
presentation of experimental data. Recently pub- 
lished data!® of benzene-methanol for 35°C are 
satisfactorily represented by means of B= 1.050; 
C=-—0.116; K=6.1, the data on carbon tetra- 
chloride-methanol by means of B=1.105; C= 
—0.148; K=6.1. The deviations of the mole 
fraction of the vapor and of the pressure are of 
the same order as those obtained by Scatchard, 
Wood and Mochel with three independent coef- 
ficients for each system. 

The characteristic influence of association is 
more or less clearly noticeable in the data 
available for benzene-phenol," benzene-aniline 
and benzene-benzonitrile,” benzene-propanol,’* 
and in various systems with ethanol and butanol." 

8 J. B. Ferguson, J. Phys. Chem. 36, 1123 (1932). 

9]. B. Ferguson, M. Freed, and A. C. Morris, J. Phys. 
Chem. 37, 87 (1933). 

10 G, Scatchard, S. E. Wood, and J. M. Mochel, J. Am. 


Chem. Soc. 68, 1957, 1960 (1946). 
“a A. R. Martin and C. M. George, J. Chem. Soc. 1933, 
14. 
12 A, R. Martin and B. Collie, J. Chem. Soc. 1932, 2662. 
18 Siang Chieh Lee, J. Phys. Chem. 35, 3558 (1931). 
4 C, P. Smyth and E. W. Engel, J. Am. Chem. Soc. 51, 
2660 (1929). 


The association indicated by the vapor pres- 
sures of the system triethylamine-water” is 
undoubtedly due to water since triethylamine is 
not appreciably associated.'® 

A qualitative check may be seen in the fact 
that the values of B are very large and that the 
values of K, with the exception of K=30 for 
hexane-methanol, are reasonable. No term similar 
to the association function A has been found in 
any system of non-associating components. It 
may be concluded, therefore, that this function 
actually is characteristic of associating systems. 

According to (10) and (19), the mole fraction 
of the monomer in the pure alcohol is 


Ny=1/(1+K). (49) 


Since this quantity as a property of the pure 
component cannot depend on the nature of the 
other component, the present model requires 
that the association constant be a property of the 
alcohol. 

Since K is extremely sensitive against experi- 
mental errors, only very accurate data can con- 
tribute to a test of this conclusion. The data of 
Benedict et al.? were obtained at atmospheric 
pressure, i.e., at varying temperatures. However, 
most of the measurements on toluene-methanol 
(K =3.8) were carried out in the neighborhood 
of 64°, those on heptane-methanol (K = 6.6) close 
to 58°. 

The association constants for methanol de- 
rived at 64° (with toluene) and 35° (with benzene 
and carbon tetrachloride) lead to the reasonable 
value of 3300 cal./mole for the heat of associ- 
ation. The value K=4.4 resulting for 55° has 
been used to represent the data obtained by 
Scatchard, Wood and Mochel'® at this tem- 
perature for the systems with benzene (B=0.981 ; 
C= —0.113) and carbon tetrachloride (B = 1.027; 

= —0.120). Here, too, the deviations of the 
mole fraction of the vapor and of the pressure 
are of the same order as in the equations of 
Scatchard e¢ al. with three independent coef- 
ficients. 

The interpolation between 35 and 64° leads 
to K=5.1 at 45° and K=4.2 at 58°. The devia- 
tions of the values 30 (with hexane) and 6.6 (with 


15L. D. Roberts and J. E. Mayer, J. Chem. Phys. 9, 
852 (1941). 
16 R. T. Lattey, J. Am. Chem. Soc. 91, 1959 (1907). 
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heptane) derived directly from experimental data 
for these two temperatures are probably due to 
errors in the measurements. 

The best test of the present results would be a 
comparison of the association constants derived 
from vapor pressures with values obtained from 
the intensity of infra-red absorption bands of the 
hydroxyl group. At present no data are available 
for this comparison. Further tests could be made ° 
in connection with data on dielectric constants 
(cf. Wilson’), the frequencies of the hydroxyl 
bands (cf. Rasmussen®) and with data on heats 
of mixing. 


7. SYSTEMS OF TWO ASSOCIATING COMPOUNDS 


If both components associate, one has, in the 
general case, to consider bonds between mole- 
cules of the first component (equilibrium con- 
stant K’), between molecules of the second com- 
ponent (K’’), from a molecule of the first kind to 
one of the second (K"’’) and conversely (K‘), 
Even with the far-reaching assumptions intro- 
duced in the first section, the general case leads 
to difficulties of such extent that a thorough 
treatment does not appear to be worth while. But 
quantitative results can be easily derived for 
two special cases. They are limiting cases of such 
nature that the whole field appears to be quali- 
tatively covered. . 

If the association constants for all four types 
of bonds are equal, K’=K”=K’’=Kiv, no 
association term at all appears in the free energy. 
This somewhat surprising result seems to be in 
satisfactory agreement with data for systems 
both components of which tend equally to the 
formation of hydrogen bonds like systems of two 
alcohols. 

The second special case is characterized by the 
absence of interassociation between the com- 
ponents, K’” = Kiv=0. One will expect that this 
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case is realized only by a kind of association 
entirely different from that caused by hydrogen 
bonds. A calculation which runs strictly parallel 
to the course followed in the third and fourth 
section of this paper furnishes again Eqs. (37) 
and (42) with 


K=K'+K"+K'K". (50) 


In other words, a system of this kind behaves 
just like a system with a single associating com- 
ponent, except that the “apparent” association 
constant has a different significance. 

The system benzene-cyclohexane appears to be 
an example of this kind. Scatchard, Wood and 
Mochel'? have already interpreted their measure- 
ments by means of association. The association 
is so small in this case that the influence of the 
volumes as indicated by (38) is still significant. 
The approximation (46) is valid. The experi- 
mental data are satisfactorily represented by 
(38) with K=0.244 (40°C) 
and K =0.184 (70°C). 

The results obtained so far indicate a fairly 
simple general interpretation of the thermo- 
dynamic properties on non-electrolyte solutions. 
The B-term in series (1) is known to represent in- 
teraction and the C-term indicates in the simplest 
cases the influence of the difference of the 
molal volumes. The D-term indicates the in- 
fluence of association and is better replaced by 
the association function A. This interpretation 
may offer a suitable starting point for a discussion 
of the coefficients of series (1) in relation with the 
nature of the components. 

The authors wish to express their thanks to 
Dr. R. W. Millar, Dr. R. S. Rasmussen and Dr. 
J. N. Wilson for valuable discussions and advice. 


17G. Scatchard, S. E. Wood, and J. M. Mochel, J. 
Phys. Chem. 43, 119 (1939). 


res- 2 
is 
e is 
fact 
the 
for 
1 in 
tion 
ms. 
tion 

49) 
ure 
the 
ires 
the 
eri- 
on- 
of 
eric 
yer, 
nol 
ood 
ose 
de- 
ene 
ble 
has 
by 
81; 
the 
ure 
of 
ef - 
ads 
ria- 4 
ith 


THE JOURNAL OF CHEMICAL PHYSICS VOLUME 15, NUMBER 12 


DECEMBER, 1947 


Foreign Ion Rejection in the Growth of Sodium Chloride ane Crystals 
from the Melt* 


RayMonD H. McFEE** 
‘George Eastman Research Laboratory of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received June 20, 1947) 


A study is made of the process of purification that takes place when single crystals are 
grown from the melt. The rejection of potassium and copper ions from solid solution in sodium 
chloride crystals is measured quantitatively by analyzing by spectrochemical means sections 
of cylindrical crystals of sodium chloride to which has been added a known quantity of im- 
purity. The theory of the distribution of impurity over the height of the crystal is derived in 
terms of the initial impurity concentration, height of completed crystal, and the fraction of 
the impurity concentration in the melt included in the crystal, which fraction is called the 
purification coefficient. Good agreement is obtained between theory and experiment for an 
initial concentration of 0.1 atomic percent potassium in sodium chloride. Purification coefficients 
for several crystals containing potassium and a few crystals containing copper are measured 
and correlated with rate of growth and temperature gradient conditions. 


INTRODUCTION 


GOOD single crystal which is to be used 

for optical purposes is characterized by its 
lack of excessive absorption of radiation in the 
optical spectrum. The increased absorption which 
is exhibited by inferior optical crystals is due, 
in general, to the presence of distributed im- 
purities,*** which may consist of ions of a 
foreign material in solid solution or of neutral 
atoms of one of the ionic constituents. The 
growth conditions for an optical crystal are 
usually adjusted experimentally to give maxi- 
mum rejection of these impurities. Stockbarger! 
found that in order to produce optical lithium 
fluoride crystals extremely pure raw material 
must be used, and the crystallizing conditions 
adjusted very carefully. It was observed that 
additional impurity material appeared at the 
top of the crystal after growth from the bottom 
upward. The phenomenon of the purification of 
materials by recrystallization from the melt has 
been utilized before. Davey? mentions that cer- 
tain crystalline materials have been purified by 
fusing and allowing to recrystallize slowly. 

* Taken from a thesis submitted in partial fulfilment of 
the requirements for the degree of Doctor of Philosophy 
- eae from the Massachusetts Institute of Technology, 

** Now with Photoswitch Incorporated, Cambridge, 
Massachusetts. 

*** Absorption caused by “F’”’ centers is not ordinarily 
observed in freshly prepared synthetic crystals. 

1D. C. Stockbarger, Rev. Sci. Inst. 7, 133 (1936). 


2W. P. Davey, Crystal Structure and its A pplications 
(McGraw-Hill Book Company, Inc., New York, 1934). 


In order to make some quantitative measure- 
ments of the rejection of foreign ions in solid 
solution in the process of the growth of an ionic 
crystal a representative crystal was selected 
which would lend itself to the method of study. 
Sodium chloride was chosen because of the ease 
in growing single crystals from C.P. grade ma- 
terial, and its water solubility, which would 
facilitate the analysis of the crystals. Two 
foreign ions were selected as impurities: potas- 
sium and bivalent copper, both in the form of 
the chloride. The selection of impurity ions was 
also restricted by consideration of the analytical 
problem, which in this case involved the spectro- 
chemical method. 


EXPERIMENTAL METHOD 
Growth of Crystals 


The general method of Bridgman,’ as modified 
by Stockbarger,* was used to grow the single 
crystals. Twenty-five grams of C.P. sodium 
chloride, to which was added an impurity 
chloride in the proportion 0.1 atomic percent 
impurity to sodium, was crystallized in a conical 
bottom platinum crucible 19 mm in diameter 
and 75 mm long. The completed crystal was 19 
mm in diameter and 44 mm long. Crystallization 
was accomplished by melting the powdered 
stock in a region of the furnace kept at a temper- 
ature above the melting point, then lowering the 


3 P, W. Bridgman, Proc. Am. Acad. 60, 305 (1925). 
4D. C. Stockbarger, J. Opt. Soc. Am. 27, 416 (1937). 
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crucible containing the melt slowly through an 
annular baffle into a region of the furnace which 
was below the melting point, causing the crystal 
to grow upward from the pointed bottom. 
Temperatures were adjusted so that crystalliza- 
tion took place close to the level of the baffle. 
The crystallizing arrangement is shown Fig. 1. 
The temperature gradient existing in the melt 
near the crystallizing level was measured by 
means of a platinum-platinum, 10 percent rho- 
dium thermocouple situated at a fixed position 
inside the crucible (Fig. 1, £). Temperature 
readings were taken at intervals, as the thermo- 
couple and crucible were lowered at a constant 
rate, until just before crystallization took place 
at the level of the thermocouple junction, when 
the thermocouple was extracted. The tempera- 
ture was found to decrease uniformly from about 
50°C above the crystallizing point, then to 
decrease uniformly at a slightly greater rate in 
the crystal. The effective temperature gradient 
was taken as the quotient of the slope of the 
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Fic. 1. Crystal growing furnace: A =thermocouple 
depth gauge; B = fire brick insulation; C=alundum furnace 
core with No. 16 Nichrome heater winding; D=platinum 
crucible; E=gradient thermocouple; F=furnace thermo- 
couples; G =baffle; H=melt; J=crystal. 
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temperature vs. time plot in the melt near the 
crystallizing point and the rate of growth. 

The rate of growth was determined by the 
rate at which the crucible assembly was lowered 
by the elevator mechanism. This mechanism 
consisted of an adjustable gear train driven by 
a small synchronous motor and connected to a 
rack on the elevator shaft. 

The level of crystallization was determined by 
means of a platinum wire probe which dipped 
into the melt through the top. With the furnace 
arrangement shown in Fig. 1 it was found that 
the position of the crystallizing level relative to 
the baffle did not change appreciably during the 
growth of the crystal when the furnace tempera- 
tures were properly adjusted. 

The furnace temperature was controlled by 
means of a special temperature regulator® oper- 
ating from the output of the thermocouple 
projecting into the upper section of the furnace. 
The temperature at this point was maintained 
constant within 0.1°C. 

The crucible containing the completed crystal 
was taken out of the furnace, inverted, and the 
crystal melted out by lowering into the hot 
furnace again, after the technique of Stock- 
barger.® 


Analysis of Crystals 


The cylindrical single crystals were placed in a 
jig which held the crystal while it was sliced 
perpendicular to the cylinder axis. Two mm 
thick slices were extracted at separations of 


5D. C. Stockbarger, Rev. Sci. Inst. 10, 205 (1939), 
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TABLE I. 

Crystal Rate of Temp. 
No. e growth grad. an an/Co A A AR 
56 0.073 1.172 os 0.0174 0.167 0.155 0.157 0.178 

0.22 0.0203 0.195 0.150 
0.378 0.0256 0.246 0.165 
0.530 0.0310 0.298 0.158 
0.694 0.0449 0.431 0.156 
0.850 0.0800 0.769 0.156 
68 0.378 1.172 3.83 0.0276 0.265 0.181 0.193 
0.530 0.0387 0.372 0.205 
69 0.378 1.172 7.32 0.0252 0.242 0.165 0.179 
0.530 0.0365 0.351 0.192 
55 0.229 2.344 — 0.0210 0.202 0.160 0.155 0.175 
0.378 0.0232 0.223 0.151 
0.530 0.0303 0.291 0.156 
0.694 0.0415 0.399 0.146 
0.850 0.0830 0.798 0.164 
57 0.073 2.344 ~ 0.0192 0.185 0.175 0.163 
0.229 0.0212 0.204 0.162 
0.378 0.0261 0.251 0.172 
0.530 0.0341 0.328 0.178 
0.694 0.0460 0.442 0.165 
0.850 0.0670 0.644 0,123 
70 0.378 2.344 4.65 0.0313 0.301 0.209 0.208 
0.530 0.0388 0.373 0.206 
54 0.073 4.688 _ 0.0203 0.195 0.185 0.173 0.196 
0.229 0.0226 0.217 0.171 
0.378 0.0260 0.250 0.170 
0.530 0.0331 0.318 0.173 
0.694 0.0489 0.470 0.179 
0.850 0.0820 0.788 0.161 
58 0.073 4.688 —- 0.0277 0.266 0.251 0.206 
0.229 0.0245 0.236 0.189 
0.378 0.0300 0.288 0.198 
0.530 0.0440 0.423 0.240 
0.694 0.0520 0.500 0.194 
0.850 0.0820 0.788 0.161 
59 0.073 4.688 — 0.0200 0.192 0.177 0.186 
0.229 0.0230 0.221 0.174 
0.378 0.0302 0.290 0.198 
0.530 0.0376 0.361 0,195 
0.694 0.0519 0.499 0,194 
0.850 0.0880 0.845 0.180 
60 0.378 4.688 -~ 0.0310 0.298 0.205 0.205 
61 0.378 4.688 _ 0.0280 0.269 0.185 0.185 
62 0.378 4.688 _— 0.0320 0.308 0.212 0.212 
65 0.378 4.688 7.68 0.0301 0.289 0.200 0.210 
0.530 0.0410 0.396 0,220 
66 0.378 4.688 3.26 0.0266 0.256 0.174 0.192 
0.530 0.0393 0.202 0.210 
64 0.378 9.375 2.56 0.0336 0.323 0.224 0.229 0.257 
0.530 0.0440 0.423 0.240 
67 0.378 9.375 7.20 0.0288 0.276 0.190 0.215 
0.530 0.0440 0.423 0.240 
72 0.378 9.375 7.72 0.0409 0.392 0.278 0.297 
0.530 0.0550 0.529 0.315 
73 «0.378 9.375 6.58 0.0360 0.346 0.242 0.268 
0.530 0.0520 0.500 0.294 
81 0.530 9.375 4.14 0.0498 0.479 0.278 0.278 
78 0.378 18.75 1.42 0.0509 0.489 0.360 0.360 0.375 
82 0.530 18.75 7.23 0.0578 0.555 0.337 0.337 
84 0.530 18.75 3.05 0.0678 0.650 0.423 0.423 
85 0.530 18.75 5.03 0.0667 0.641 0.411 0.411 
86 0.530 18.75 4.64 0.0593 0.570 0.344 0.344 
87 0.530 25.00 5.37. 0.0920 0.885 0.715 0.715 0.715 


nN 
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about six mm. The concentration of impurity in 
a slice was taken to represent the concentration 
at the average height of the slice in the crystal. 

Spectrochemical means were used to determine 
the impurity concentration in the crystal slices. 
The concentrations were in such a range as to be 
favorable to this type of microanalysis. 

Analyses for potassium in the sodium chloride 
crystals were made by photographing the spec- 
trum of the material excited by a d.c. are through 
a logarithmic step-sector disk. Rubidium chloride 
was added to the salt as an internal standard for 
the analysis. Microphotometer measurements of 
sensitive lines of potassium at \7665A and 
rubidium at \7800A recorded on infrared sensi- 
tive film were used as indications of the concen- 
trations.°® 

Concentrations of copper ions were measured 
by exciting the samples by means of an a.c. high 
voltage arc. Silver was used as the internal 
standard in this case. The copper line at \3274A 
and the silver Jine at \3281A were used for the 
analyses. 


THEORY 
Distribution of Impurity Ions in the Crystal 


Crystallizing conditions in the above described 
apparatus are such that the following assump- 
tions may be made, particularly for the single 
crystals grown at low rates: (1) the rate of 
growth and temperature gradient remain con- 
stant throughout the growth; (2) the impurity 
ions are distributed uniformly throughout the 
melt; (3) the foreign ions are contained in solid 


6L. W. Strock, J. Opt. Soc. Am. 32, 103 (1942). 
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solution in the crystal; (4) the phase boundary 
between melt and crystal is a plane perpendicular 
to the cylinder axis. Figure 2a represents the 
picture of the crystal growth. The height dimen- 
sion is measured from a level one third the 
altitude of the conical bottom section below the 
shoulder, corresponding to the bottom of a 
cylindrical section of equivalent volume. 

With constant crystallizing conditions, the 
probability that an impurity ion will become 
trapped in the crystal lattice is constant. If the 
concentration of impurity ions in the melt is 
small, the fraction of the number of impurity 
ions per unit mass in the melt that is included in 
the crystal is a constant. Then the concentration 
of impurity ions in the crystal, a,, at a height h 
is proportional to the concentration of impurity 
ions in the melt, b,, at the point when the crystal 
has grown to a height h, or 


a,=Abd,, (1) 


where A, the constant of proportionality, may 
be called the purification coefficient. This 
quantity never exceeds unity. 

From this relation the distribution of impurity 
ions in a crystal grown under the above crystal- 
lizing conditions may be derived. Let dN, be 
the change in the total number of impurity ions 
in the melt when the crystal grows from a height 
h to a height h+dh. The number of impurity 
ions per unit height in the crystal at height h is 
then —dN,,/dh. lf a, and b, are expressed in 
atomic percent impurity ions to sodium ions, 


1007 dN» 
a,=—( -—}, 

100 7 Nn 
n-—(—), (3) 

\ 


where , is the number of sodium ions per unit 
height in the crystal, 2» is the number of sodium 
ions per unit height in the melt, NV, is the total 
number of impurity ions in the melt, and h,, is 
the height of the melt. 

hm bears the relation to H,, the height of the 
completed crystal, 


Ne 


Nm 


TABLE II. 
Crystal Rate of Temp. : = 
No. h/He growth grad. ah an/Co A A 


94 0.530 4.688 1.71 0.0315 0.315 0.170 0.170 
75 0.073 4.688 6.40 0.0107 0.107 0.095 0.0644 


0.229 0.0094 0.094 0.070 
0.378 0.0100 0.100 0.062 
0.530 0.0145 0.145 0.069 
0.694 0.0191 0.191 0.061 
0.850 0.0365 0.365 0.060 


93 0.530 9.375 


61 0.0518 0.518 0.307 0.307 
91 0.530 9.375 7 


‘70 «0.305 


aN 


By substitution and integration 


5 
H, ©) 


100 h 

where JN is the total number of impurity ions 
present, or 


’ (6) 


c 


where Cy is the concentration of impurity ions 
initially present in atomic percent impurity to 
sodium ions. 

Figure 3 shows a plot of the fractional im- 
purity concentration, a,/Co, as a function of A 
for a series of values of h/H,. corresponding to 
those values for the six sample slices analyzed 
(Fig. 2b). If a cross section of these curves is 
taken at a given value of A, the resulting curve 
represents the relative distribution of impurities 
throughout the height of a crystal grown under 
conditions such that the purification coefficient 
is A. 


EXPERIMENTAL RESULTS 


The results of measurements taken on twenty- 
five sodium chloride crystals containing initially 


| || 
CRYSTAL NO 56, KIN NeCi| | | 
= | @ CRYSTAL NO. 59, K IN NaC 
} 1@ CRYSTAL NO 7S, Cu IN NaC! | 
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0.104 atomic percent potassium (0.100 atomic 
percent added plus 0.004 atomic percent potas- 
sium impurity) are given in Table I. Each 
crystal used for analysis was homogeneous, single, 
clear, and colorless. Each measurement of the 
potassium concentration, @,, represents the aver- 
age of five spectrographic determinations. Tem- 
perature gradient measurements were taken on 
crystals from No. 64 on. A values corresponding 
to the a, figures were obtained by interpolation 
in Fig. 3 on the curve representing the level of 
the sample analyzed. a;/Cy) represents the frac- 
tional impurity concentration, that is, the frac- 
tion of the original impurity concentration 
remaining in the crystal at height 4. A is the 
average value of A for the various samples from 
a single crystal. Since no systematic variation of 
A with temperature gradient is evident, the 
average A for all crystals grown at the same 
rate, Ar, is taken. 

By way of comparison with the results ob- 
tained with potassium impurity four crystals 
containing initially 0.1 atomic percent copper 
were grown and analyzed (Table II). Crystal 
number 75 was analyzed over the complete 
height. 

The fractional impurity concentration as re- 
lated to the height in two crystals containing 
potassium and in one crystal containing copper 
are shown in Fig. 4. The points are taken from 
the experimental data as given in Tables I and 
II. The solid curves are plotted from the equation 
for a,/Co (Eq. 6), where A is taken as the 
average, A, of the A values of all six samples. 

The variation of the average purification 
coefficient, Ar, with rate of growth for the 
potassium-containing crystals, is shown in Fig. 5. 


Since the point at 25 mm per hour represents 
only one crystal, its position is not as well 
established as those of the other points. 


DISCUSSION 


It would appear from the agreement of the 
experimental points for the two potassium 
crystals with the theoretical curves in Fig. 4 
that growth conditions approached the ideal ones 
as postulated in the theory in those instances. 
The rising tendency for the points at the lower 
end of the curve taken for the copper-containing 
crystal may be due to a slight settling out of the 
heavier copper ions in the melt. As might be 
expected, the copper impurity is rejected more 
completely from the crystal than in the potas- 
sium, since the latter possesses ionic character- 
istics more similar to those of sodium. Even in 
the case of potassium impurity the foreign ion 
concentration is less than that of the raw 
material up to a thin layer at the top of the 
crystal. 

Figure 5 and Table I indicate the relation 
between the degree of~ purification and the 
crystallization conditions. Figure 4 is significant 
in that it shows that a minimum value of 
potassium concentration is reached at a finite 
rate of growth, in this case near four millimeters 
per hour. The apparent random function of the 
temperature gradient in eliminating potassium 
from sodium chloride may be due to the insuffi- 
ciently extensive range of temperature gradients 
used. However, in the case of copper impurity 
the higher gradient appears to be more effective 
in the elimination of the foreign ions, while the 
lower rate of growth gives better purification as 
in the case of potassium. Because of the greater 
dissimilarity of copper and sodium ions, the 
inclusion of the copper in the crystal is more 
critically dependent upon growth conditions. 

These data indicate that this general method 
might be applied to the preparation of higher 
purity raw materials for the growth of single 
optical crystals. This process is particularly 
adaptable for this purpose since the same equip- 
ment and techniques may be used for the pre- 
liminary purification of the material as is used 
for the growth of the crystal itself. Many sub- 
stances used for optical crystals are quite difficult 
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to prepare of sufficient purity by chemical means, ACKNOWLEDGMENT 


notably lithium fluoride! and calcium fluoride. The writer would like to express his apprecia- 
This method offers a means of preparing a fairly tion for the advice and encouragement of Pro- 
large amount of purer material in a relatively fessor D. C. Stockbarger, who suggested the 
short time. problem. 
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Theory of Burning Velocity. II. The Square Root Law for Burning Velocity! 


CHARLES TANFORD? AND RoBert N. PEASE 
Department of Chemistry, Princeton University, Princeton, New Jersey 


(Received June 27, 1947) 


Previous papers by the authors having suggested that the diffusion of active particles 
(chiefly hydrogen atoms) from the flame front is the controlling factor in combustion in 
Bunsen-type burners, an. equation for burning velocity based on such a concept is derived. 
This equation shows that the burning velocity should be proportional to the square root of an 
expression in which the most important factor is the sum of the products p;D; for all effective 
atoms or free radicals, ; being the concentration of an atom or radical at the flame front, and 
D; its coefficient of diffusion into cold, unburnt gas. The burning velocity equation is applied 
to the combustion of carbon monoxide and hydrogen, and excellent agreement between cal- 
culated and experimental values of the burning velocity is obtained. The equation is also used 
to account for the effect of pressure upon burning velocity. 


INTRODUCTION 


REVIOUS papers by the authors on Bunsen- 

type burner flames (with laminar flow)*4 
have shown that there exists in moist carbon 
monoxide flames a close correlation between the 
calculated equilibrium hydrogen atom concen- 
trations at the flame front and the experimentally 
determined burning velocity, and have suggested 
that this correlation can be explained if it is 
assumed that the diffusion of active particles 
(chiefly hydrogen atomis) from the flame front 
into the unburnt gas controls the combustion 
process. In this paper an equation for burning 
velocity based on this idea is derived, and is 
given a number of practical applications. 


THE SQUARE ROOT LAW FOR 
BURNING VELOCITY 


As in our Paper I on burning velocity,’ we 
consider the simple model of the combustion 


1 Abstracted in part from a thesis submitted by Charles 
Tanford in partial fulfilment of the requirements for the 
degree of Doctor of Philosophy at Princeton University. 

? National Research Council Predoctoral Fellow in 
Chemistry, 1946-47. 

(94s) Tanford and R. N. Pease, J. Chem. Phys. 15, 431 

4C. Tanford, J. Chem. Phys. 15, 433 (1947). 


zone illustrated by Fig. 1. The flame front is 
defined as the surface where combustion has 
proceeded to equilibrium, and is used as the 
origin of a fixed coordinate system. The unburnt 
gas moves in the direction of the arrow, the 
flow being taken normal to the flame front, so 
that in the cold unburnt gas the flow velocity 
must be the burning velocity as defined in the 
usual manner. In order to arrive at soluble 
differential equations, it is necessary to assume 
that the whole combustion zone is at one con- 
stant mean temperature, 7;,, which is, in general, 
taken to be some fraction of the flame tempera- 
ture, 7;. The effect of this simplification is to 
give to the density and the flow velocity values 
Pm and um, which are constant over the entire 
combustion zone. It should be pointed out that 
the actual value chosen for the mean tempera- 
ture can affect the predicted burning velocity 
only to a small extent. 

We now consider a portion of the flame front 
of area dS and a cylinder of this same cross- 
sectional area going from it into the unburnt gas. 
It is assumed that ignition of unburnt gas can 
take place at any point in this cylinder by 
means of active particles (chiefly hydrogen 
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FLAME 
FRONT 
FLOW OF 
UNBURNT GAS 
ZleREACTION 1ONE> 


*=Q —DISTANCE, x ———> 


Fic. 1. Model of flame front and combustion zone. 


atoms) present there as a result of diffusion from 
the flame front. The concentration in molecules 
per cc of any such active particle at any point 
has already been obtained in Paper I, and is 
given by 


Bum 
‘ Ci=C1; exp(- (1) 


mi 


where ¢;; is the equilibrium concentration of the 
active particle at the flame front, Dn; is the 
diffusion coefficient of the active particle into 
unburnt gas at the mean temperature 7, and 
B; is a factor, usually close to unity, which 
corrects for the loss of radicals due to chemical 
processes. The evaluation of this last factor for 
hydrogen atoms is discussed in detail in Paper I. 

Now, the rate of formation of combustion 
product at any point must certainly depend upon 
this active particle concentration. It can be 
expected to be the sum of a number of terms, 
one for each effective atom or free radical. The 
concentrations of species other than these active 
particles can also be expected to enter into the 
expression. For the sake of simplicity we shall 
assume here that of the latter only the concen- 
tration of combustible, C, is of importance, and 
that the rate of formation of combustion product 
in molecules per cc per sec. at any point can be 
represented by 


d(product) /dt =>): kiciC 


where k; is a rate constant appropriate to the ith 
active species. Then the total rate of formation 
of product molecules in the cylinder of Fig. 1 in 
molecules/sec. is 


f kic;CdSdx. 
r=0 


Now, from the initial composition of the 
combustion mixture it is possible to compute the 
total rate of formation of combustion product 


PEASE 


which is required to bring about complete com- 
bustion, and this rate must be equated to the 
above expression. To calculate this rate we 
consider the volume of gas (measured at 298°K) 
which crosses the section dS of the flame front 
per second. This volume is given by the product 
of the burning velocity (uo) and the cross- 
sectional area (dS). If Q is the mole fraction of 
potential combustion product in this gas,® then, 
where L is the number of molecules in one cc of 
gas at 298°K, the number of molecules of product 
required to be formed per second over the cross- 
sectional area dS for complete combustion must 
be LQuodS. Since, by definition, combustion 
must be complete at the flame front, we then 
have 


z=0 


We now substitute for ¢; from Eq. 1, and 
assume a mean value independent of x for the 
concentration of combustible, so that we obtain 


Bittn 
LQ ud S= ae 
i 0 


mi 


=2, (2) 


Um 


The equilibrium composition of the burnt mix- 
ture is usually in terms of partial pressures or 
mole fractions,® so that we replace ¢1; by Lpi/On 
where 6, is T/T, the ratio of the mean temper- 
ature of the combustion zone to the temperature 
of the unburnt gas (usually 298°K), and L/@» 
therefore the number of molecules in one cc of 
gas at the mean temperature. Furthermore, the 
diffusion coefficient varies approximately as the 
square of the temperature,® so that we can 
replace Dn; by Di§n?, where D; is the diffusion 
coefficient of the ith active species into cold, 


5 The compositions of the unburnt gas and of the final 
equilibrium mixture at the flame front are most conveni- 
ently expressed as volumes percent, or, better as mole 
fractions. We shall use the symbols g, to represent mole 
fractions in the unburnt gas, and the symbols p; to 
represent calculated equilibrium mole fractions at the 
flame front. If the flame is at a pressure of one atmosphere, 
these mole fractions are, of course, equivalent to partial 
pressures. 

6S. Chapman and T. G. Cowling, The Mathematical 
Theory of Non- — Gases (Cambridge University 
Press, 1939), p. 248. 
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unburnt gas. Finally, the flow velocity, uw», can 
be related to the burning velocity, i.e. the flow 
velocity in cold, unburnt gas, by the relation 
Um = U9 m. 

Making these substitutions in Eq. 2, we find 
that L, @,, and dS drop out, and we obtain the 
square root law for burning velocity, 


(3) 
i OB; 


Since the existence of combustion limits im- 
poses a relatively narrow range of variability 
upon the concentration of combustible, C, and 
upon the mole fraction of potential combustion 
product, Q, and since the factors B; vary but 
little from mixture to mixture, we can say that, 
for the purpose of comparing the burning ve- 
locities of different mixtures in which the same 
kinetics (i.e., the same k,;) are likely to be 
involved, the burning velocity must be roughly 
proportional to (>>; p,)*. This immediately 
explains the correlation between burning velocity 
and hydrogen atom concentration found for 
mixtures containing carbon monoxide.* 


APPLICATION TO CARBON MONOXIDE 
COMBUSTION 


In the combustion of carbon monoxide the 
mole fraction of potential combustion product in 
the initial mixture is the mole fraction of carbon 
monoxide, or (since one mole of oxygen can 
oxidize two moles of carbon monoxide) twice the 
mole fraction of oxygen, whichever of these is 
smaller. This is conveniently represented by the 
symbol (gco, 2g02) <. If the average mole fraction 
of carbon monoxide over the combustion zone 
is, say, 0.7 of that present initially, the average 
concentration of combustible, C, can be ex- 
pressed in terms of that initial concentration. 
For, since the number of molecules in one cc of 
gas at the temperature T,, is L/@m, we then have 


C=0. 7L4co/Om.- 


Considerable evidence exists which suggests that 
oxygen atoms are not of importance in initiating 
carbon monoxide combustion (for example, com- 
bustion in the absence of moisture, even with an 
excess of oxygen, proceeds only with the greatest 
difficulty), so that it appears reasonable to sup- 


pose that only hydrogen atoms and hydroxyl 
radicals of the available active particles are of 
interest to us. Of these, the contribution of 
hydroxyl is usually small, because Don is small. 
In a few of the mixtures for which calculations 
have been made, however, the concentration of 
hydroxyl is very much larger than that of 
hydrogen atoms, and in those cases good agree- 
ment of experimental burning velocities with 
those predicted by Eq. 3 can be obtained only 
by inclusion of the part played by hydroxyl. 
Thus the special form of this equation for carbon 
monoxide combustion becomes 


0.7L dco 
Uy = 
(dco 29q02)< 


kupuDyu koubonDon\ 
x( + ) . (4) 
By Bou 


This equation can be simplified to some extent. 
The correction term Box is found to be always 
within one or two percent of unity, and can 
therefore be neglected. Further, in all mixtures 
containing chiefly carbon monoxide, oxygen and 
nitrogen, in any proportion whatsoever, diffusion 
coefficients have about the same value. At 298°K 
Dy is 1.95 cm?/sec. (see paper I) and Don is 
0.29 cm?/sec. Finally, we can assume that the 
rate constant k has about the same value for H 
as for OH. Placing the mean temperature, 7°, 
equal to, say, 0.7 of the equilibrium flame 
temperature, and substituting for Z its numerical 
value, 2.48 10!® (298°K, one atmosphere pres- 


TABLE I. Burning velocities in carbon monoxide flames. 


1 Py Pou uo, 
Ti X104 cm/sec. 


Composition “7, atmos. atmos. By calc. obs.° 


35%CO,® 8.5%O2, 56.5%N2 5.8 O01 — 1.08 22 25 
26.5%CO,* 44.1%02,294%N2 78 0.7 206 1.23 36 30 
12.6% 02, 474%N2 78 25 34 132 37 30 
40%CO.* 474% N2 7.8 3.8 49 1.28 41 40 
59.1%02, 0.9%N2 91 56 56 125 Si 48 
40%CO,* 24%02, 36%N2 8.9 116 56 139 56 St 
60% CO,» 39.4%02, 0.6%N2 98 24 86 135 67 68 
60%CO* 24%0s, 16%N: 96 38 90 139 81 83 
60%CO,* 39.4%02, 0.6%N: 98 39 134 1.32 81 93 
70%CO.* 24%02, 6%N: 9.7 52 102 137 96 98 
80%CO,* 19.7%02, 0.3%N2 96 59 78 1.37 112 106 


® CO contains 1.35% H2O and 1.5% Ho. 
> CO contains 1.35% H20O. 
as = Jahn, Der Ziindvorgang in Gasgemischen (Oldenbourg, Berlin, 
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TABLE II. Burning velocities in hydrogen flames. 


Pu Pou 0, 
X104 cm/sec 

Composition T) atmos. atmos. Bu calc, obs.* 
40%H2, 5%O2, 52.5%N2 5.0 0.1 1.002 120 100 
8.4%O2, 31.6%Ne $3 0.6 1.04 140 175 
40% He, 10. 49.5%Ne2 6.5 5.6 04 1.16 180 200 
40%He 15%0O2, 45%Ne 8.3 85 40 1.10 330 300 
12%02, 28%Ne 32 2.55 1.15 320 350 
40%H2, 24%O2, 36%Ne2 9.2 135 380 1.01 390 360 
40%H2, 36%0 24%N2 9.2 96 520 1.06 390 380 
60%H2e 16%02, 24%N 8.7 214 50 1.15 530 485 
60%H:z 24%02, 16%Ne 99 590 400 1.02 700 
60%He 32%O2, 8%Ne2z 10.1 590 890 1.02 690 730 


wear Jahn, Der Ziindvorgang in Gasgemischen (Oldenbourg, Berlin, 


sure), we then obtain 
10% 
Up = 
(qco2q02)< 


(5) 


Before this equation is applied, it must be 
remembered that the burning velocity of carbon 
monoxide mixtures does not tend to zero as the 
free radical concentration becomes zero,’ i.e., 
there appears to be a small, probably fairly 
constant contribution to the burning velocity 
which is not dependent upon the free radical 
concentration. A plot of burning velocity against 
the square root of the free radical concentration 
shows that this radical-independent contribution 
is about 17 cm/sec. This value agrees closely 
with burning velocities observed by Ubbelohde 
and Dommer® for carbon monoxide—air mixtures 
in the complete absence of moisture or hydrogen, 
i.e., under conditions where H and OH could not 
be formed. This contribution to the burning 
velocity must be added to Eq. 5, so that we 
obtain 


4.8410" dco 


uyp=17+ 
6; (qco:2q02)< 


Pu 
x(+0.15p0n)| . (6) 
Br 


Equation 6 has been used to compute burning 
velocities of a number of mixtures containing 


3, Fig. 1. 
9 L. Ubbelohde O. Dommer, J. Gasbel. 57, 757 
1914 
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carbon monoxide. The constant & is, of course, 
arbitrary, best fit of the experimental data being 
obtained with the value 1.64x10-"%. Methods 
for the calculation of the free radical concentra- 
tions at the flame front and of the term By 
have been described in previous papers. The 
results of the calculations are shown in Table |, 
and it is seen that excellent agreement exists 
between the calculated burning velocities and 
experimental values determined by Jahn.® It 
was found that somewhat less good agreement 
is obtained for extremely rich mixtures (ratio of 
CO and O, greater than 4), and such mixtures 
have not been included in Table |. 

It can be seen at once that the value found for 
k is of the right order of magnitude. For, con- 
verting from the units of cc molecule sec. to 
the more familiar units of cc mole sec.-! by 
multiplication by Avogadro’s numer, we obtain 
k=1.00X10" cc mole sec.'. This value is 
entirely reasonable for the high temperatures 
obtaining in the combustion zone,'® and suggests 
that the activation energy of the rate-determin- 
ing step in the combustion mechanism is about 
7 or 8 kcal./mole. The fact that the same value 
for the rate constant can be used throughout, 
despite the fact that there is variation in the 
mean temperature, also requires that the activa- 
tion be of this order or lower. 


APPLICATION TO THE COMBUSTION OF 
HYDROGEN 


Assuming that the mechanism of hydrogen 
combustion is similar to that of carbon monoxide 
combustion, we obtain for the burning velocity 
an equation similar to Eq. 4. Again it is found 
that there is a contribution independent of the 
concentration of active particles, this time of 
about 110 cm/sec., so that the burning velocity 
equation becomes 


2.48 X10"kDy 


Uy = 110+| 
0; (que: 2q02)< 


pu 
—+0.21p0 
x(F+0 21p «)} (7) 


9G. Jahn, Der Ziindvorgang in Gasgemischen (Olden- 
bourg, Berlin, 1934). 

10 For reaction at every collision and zero entropy of 
activation, the bimolecular rate constant is about 10" cc 
sec.—1, 
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It should be noted that while the ratio Don/Da 
remains about constant, the value of Dy now 
depends upon the composition, being 2.35 for 
mixtures containing 40 percent hydrogen, and 
2.6 for mixtures containing 60 percent hydrogen. 

Equation 7 has been used to compute the 
burning velocities of a number of mixtures con- 
taining hydrogen, with the best value of k now 
6.5 X10-" (corresponding to a value of 3.910" 
cc mole sec.~!), i.e., slightly higher than for 
carbon monoxide. The computations are shown 
in Table II, and again good agreement is found 
between calculated burning velocities and those 
observed experimentally by Jahn.® 


THE EFFECT OF PRESSURE ON BURNING 
VELOCITY 


Equation 3 should prove to be very useful in 
predicting the effect of a number of variables, 
such as pressure, preheating, nature of diluent 
gas, etc., upon burning velocity. As an example 
we shall give here a semi-quantitative estimation 
of the effect of pressure upon burning velocities 
in general. 

We observe at once that only three of the 
factors of Eq. 3 depend upon pressure, these 
being C, p;, and D;. Moreover, for a mixture of 
given initial percentage composition, the concen- 
tration of combustible, C, in molecules per cc, is 
directly proportional to the pressure, while the 
diffusion coefficients, D;, are inversely propor- 
tional to the pressure. These two effects therefore 
cancel, and it remains to consider only the 
pressure-dependence of the free radical mole 


fractions, p;. Since these free radicals are products 
of dissociation, their mole fractions will decrease 
with increasing pressure, in general roughly as 
the square root of the pressure. (Slight changes 
in flame temperature and molecular concentra- 
tions always occur, so that an analytical expres- 
sion is not feasible.) Since the burning velocity 
itself is proportional to the square root of the /;, 
therefore, it should vary inversely, roughly as 
the fourth root of the pressure. Such an effect 
has indeed been observed for a number of 
hydrocarbon flames, as well as for carbon 
monoxide.” 


CONCLUSION 


It should perhaps be pointed out in conclusion 
that the very simple treatment of burning 
velocity given in this paper is successful only 
because the effects of the motion of the gas 
stream and of diffusion upon local free radical 
concentrations are relatively much greater than 
the effects of chemical processes. In other words, 
the concentration of active particles at any 
point depends very little upon the chemical 
reactions occurring in the gas. Of the latter, we 
have taken into account only the recombination 
of hydrogen atoms and the formation of HO», 
and it is very likely that the accuracy of the burn- 
ing velocity equation could be improved if full 
consideration could be given to chain branching 
and other chemical processes which in this simple 
treatment have been ignored. 

"1 L. Ubbelohde and coworkers, J. Gasbel. 59, 49 (1916), 


60, 225 (1917); L. Khitrin, Tech. Phys. U.S.S.R. 3, 926 
(1936). 
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Heat Capacities of Liquids and Vapors 


SIDNEY W. 


BENSON 


Chemical Laboratories, University of Southern California, Los Angeles, California 
(Received July 14, 1947) 


By the use of an equation previously derived for the energy of vaporization of a liquid: 


AE=k(Di" 


—D,") (in which k is a constant; n=5/3, and D; and D, are the molar densities of 


liquid and vapor, respectively), an expression is derived relating (@AE/dT)sat to the difference 
in molar heat capacities at constant pressure of the liquid and vapor (AC,). At low pressures 
this relation can be approximated by the equation AC, = (@AE/dT),+R (R is the gas constant, 
1.986 cal./mole-°A). Under these conditions InD,/dT),, and the 
resulting equation for AC, is AC,=(5/3)AE(0 InD,/dT),+1.99. From known properties of 
liquids AC, can be calculated using this last equation. These calculated values are found to be 
in good agreement with observed experimental values. 


N a recent paper! it was shown that a cage 

model of a liquid could be used to derive a 
simple analytic expression for the configurational 
energy of the liquid. If the same expression is 
applied to the configurational energy of the vapor 
in equilibrium with the liquid, the energy of 
vaporization is given by the difference between 
these configurational energies: 


AE vay = E-(liq.) — E-(gas). (1) 


Since the configurational energies obtained 
from the cage model are explicit functions only 
of the density of the particular phase involved, 
the energy of vaporization becomes an explicit 
function of these parameters only, and is given 
by the following equation : 


AE vap =k(Di"—D,"). (2) 


D, and D, are the molar densities of the liquid 
and gas phases, respectively, k is a constant for 
any given liquid, and m is a number which may 
be taken as 5/3 for all liquids.' (In the neighbor- 
hood of the critical region, the value n=2 gives 
a better fit with experimental data.) 

If we differentiate this expression for the 
energy of vaporization with respect to tempera- 
ture at saturation conditions, we will obtain the 
derivative (OAE/dT),,* which may be related 
thermodynamically to AC,, the difference in the 
heat capacities of the liquid and its saturated 
vapor. This latter equation may then be solved 
explicitly for AC, in terms of experimentally 
available properties of the liquid and vapor. 


1S. W. Benson, J. Chem. Phys. 15, 367 (1947). 
* The subscript s will refer to saturation. 


RELATION BETWEEN (0AE/0T), AND AC, 


From the first and second laws it is possible 
to derive the following equation for AC,? 


PAH 
+—(—). @ 
TAV\ OP 


Regrouping the terms in Eq. (3) we find: 


AHy InAV 
Ttr\ 


0 InAV 
+( ) 
0 InP T oT 

Equation 4 is an exact relation between AC,, 
(@AE/dT),, and the known properties of a liquid 
and its vapor. The first term on the right-hand 
side of the equation is the dominant term and 
negative, a fact in agreement with the observed 
property that the heat capacity of a liquid is 
greater than the heat capacity of its vapor. 

The third term in Eq. 4 is positive, and 
the second term is made of two partial deriva- 
tives, the first of which is positive and the second 
negative. If equations of state are available, the 
second and third terms may be explicitly evalu- 
ated and AC, calculated exactly. However, at 
low vapor pressures the vapor is very dilute and 


2F,. H. MacDougall, Thermodynamics and Chemistry, 
(John Wiley & Sons, Inc., New York, 1939) third edition, 
p. 126, contain a derivation from which Eq. 3 may be 
constructed. 
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may be assumed to obey the ideal gas laws. In 
such regions the following approximations may 
be made: 


InAV dlnV, 
/, 7, 
0d InAV dInV, 
ae 
Jy dInP 


Substituting these values into Eq. 4 we have, 


finally, 
+R (8) 
oT /, 


EVALUATION OF (0AE/0T), 


and 


~R/P. (7) 


If we take the derivative with respect to 
temperature of AE (at saturation) as given by 
Eq. 2, we have 


oD, aD, 
(—) = ) (9) 
oT oT 8 


Rearranging the terms this becomes 
(—) 
ar], D; \Nar/, 
(10) 
oe 


At the pressures which are under consideration 
(in the neighborhood of the boiling point or 
lower), (0D,/0T), is never greater than 1/10 of 
(dD,/dT),. In this same region D, is usually less 
than 1/200 of D;. We can, therefore, make the 
following approximations in this pressure region : 


AE=kD,". (12) 


and 


Making the substitutions from Eqs. 11 and 12 
into Eq. 10 we find: 


OAE 0 InD, 
(— =) (— (13) 


TABLE I. AC, of liquids and their saturated vapors. 


ACp(obs.) ACp(calc.) 


Mole t cal./mole_cal./mole 

Compound Wt. (°C) —°C —°C 
Ethyl ether 74 35 14.0 14.2 
Carbon tetrachloride 154 20 11.2 12.1 
Benzene 78 80 14.9 14.4 
n-hexane 86 80 16.5 15.6 
Sulfur dioxide 64 15 11.4 12.7 
Water 18 100 9.31 9,21 


Further substitution from Eq. 13 into 8 finally 
yield 
AC, =nAE(d InD,/dT),.+R. (14) 


COMPARISON OF EQ. 14 WITH EXPERI- 
MENTAL DATA 


Equation 14 gives us a method of calculating 
AC, at low pressures from the energy of vapor- 
ization of a liquid and its coefficient of expansion. 
If data are lacking for the energy of vaporization, 
AE may be approximated at the boiling point by 
Trouton’s rule or by Hildebrand’s rule. In any 
case, all of the quantities in the first term on the 
right-hand side of ‘Eq. 14 are either well known 
or easily measured experimentally. On the other 
hand, it is very difficult to obtain reliable 
experimental values for AC,. Heat capacities for 
liquids have been measured, but heat capacities 
at constant pressure for saturated vapors are 


‘difficult to obtain. 


In Table I experimental values of AC, are 
compared with values calculated with the aid of 
Eq. 14 (1=5/3). 

The experimental values for the heat capacities 
of liquids and vapors were taken from values 
listed in Landolt-Bornstein and the Jnternational 
Critical Tables. For the calculated values of AC,, 
Eq. 14 was used, with the value of » taken as 
5/3. Values of the energy of vaporization and 
the thermal coefficient of expansion of the liquid 
were calculated from the data of Young* and 
tables found in Landolt-Bornstein. The agree- 
ment between the observed and experimental 
values of AC, (Table I) are quite satisfactory 
and probably within the limits of error of the 
experimental values for the molar heat capacities. 
The agreement lends further support to the 
value of »=5/3 which was proposed in the 
previous paper.! 


3 Young, Proc. Roy. Dublin Soc. 12, 374 (1910). 
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The Experimental Determination of the Intensities of Infra-Red Absorption Bands. 
Ill. Carbon Dioxide, Methane, and Ethane* 


A. M. THORNDIKE** 
Mallinckrodt Chemical Laboratory, Harvard University, Cambridge, Massachusetts 


The absolute intensities of the chief infra-red absorption bands of carbon dioxide, methane, 
and ethane have been determined experimentally to be (in cycles per centimeter at N.T.P.): 


Carbon dioxide Methane Ethane 


670 cm™ band: 56010" 1306 cm band: 45010” 820 band: 73X10" 
2350 band: 860010" 3020 band: 90010" 1450 cm bands: 23010" 
2950 bands: 2200 10" 


These results are in satisfactory agreement with available data on infra-red dispersion and 
atomic polarization. They may be interpreted in terms of dipole moments, yu’, of bonds between 
vibrating atoms, and their rates of change with internuclear distance, du’ /dr. When this is done n 
for the CH bond, however, unexpected variations in these quantities are found for the different 
vibrations and molecules studied. The average values are about »’=0.4X10~%, du’/dr= 


+0.6X10~". For the CO bond the values are considerably larger. 


t 
C 
INTRODUCTION mining A; from absorption measurements made ¥ 
HE light transmitted a distance, L, through with a spectrograph of low resolving power,' and 
a gas at pressure, p, is given by Eq the application of this method to ethylene and 
(1). Here Jo(v) is the incident intensity, nitrous oxide has been discussed. The results of 
further measurements using this method will be ol 
I(v) =Io(v) expl —a(v) pL ], (1) reported in the following sections. The experi- Ir 
and a(v) is the absorption coefficient. ‘The in- ental procedure involves measuring the ap- = 
tensity of absorption is defined by Eq. (2). Here Parent value of the intensity—called the “‘inte- (I 
grated apparent absorption coefficient’’—as a 1S 
= function of the partial pressure of the absorbing m 
a(v)dv, (2) gas, with an inert (non-absorbing) gas added to a on 
sate constant total pressure. The limit of this ‘‘inte- in 
A; is the absolute intensity of the ith vibrational grated apparent absorption coefficient”’ as the ab- pé 
absorption band, the integral being taken over sorbing gas pressure approaches zerois the true in- 
the band in question. tensity.! The inert gas present serves to broaden al 
Direct measurement of the intensity is difficult | the individual lines making up the absorption band 
because a(v) and J(v) may fluctuate considerably and reduce the fluctuations in a(v) which cause 
even over the small frequency range transmitted _ the ‘integrated apparent absorption coefficient” 
by an infra-red spectrograph at any setting. A to deviate from the true intensity. 
spectrograph with a finite effective slit width The absorption of light by these infra-red 
averages J(v) over that frequency range, and it bands is due to changes in the dipole moment of 
is not possible to obtain a(v) directly from such the vibrating molecule when the atoms are dis- 
a measurement. In previous papers an extra- placed from their equilibrium positions. ‘The 
polation procedure has been described for deter- intensities are proportional to |dy/dQ|?, where 
u is the dipole moment of the molecule, and Q 
* The material in this paper is taken from a thesis by : : : a 
A. M. Thorndike, submitted to the Faculty of Arts and the normal coordinate for the vibrational mode 
Sciences of Harvard University for the degree of Doctor in question.2 The complete relationship is given 
of Philosophy in Chemical Physics in May, 1947. The in Eq. (3) 
research reported in this document was facilitated by ~ , Fic 
support extended Harvard University by the Navy ———— 
Department (Office of Naval Research) under Office of 1E. Bright Wilson, Jr. and A. J. Wells, J. Chem. Phys. 
Naval Research contract N5ori-76. 14, 578 (1946). an 
** National Research Council Predoctoral Fellow, now 2A. M. Thorndike, A. J. Wells, and E. Bright Wilson, $ 
at Brookhaven National Laboratory, Upton, New York. Jr., J. Chem. Phys. 15, 157 (1947). 46, 
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A;=(Nw/3c) { | 
+ | Oty /IQ;|?+ | (3) 


The intensities may be interpreted in terms of 
the structure of the molecule by assuming that 
the dipole moment, u, of the molecule can be 
regarded as the sum of dipole moments, u,’, of 
the bonds forming the molecule. With such a 
picture the intensities of bending modes deter- 
mine yp’, stretching modes dy’/dr.2 These quan- 
tities might reasonably be expected to be charac- 
teristic properties of any particular bond, con- 
stant for different vibrations and different 
molecules, but the results show that this is true 
only as a very rough approximation. 

In the following sections values are given for 
the intensities of the fundamental bands of COs, 
CH,, and C,Hs. These are interpreted in terms 
of the properties of CO and CH bonds. 


CARBON DIOXIDE 


The intensities of the carbon dioxide bands are 
of special interest because: (a) measurements of 
infra-red dispersion are available* which may be 
compared with our values for the intensities, and 
(b) absorption of heat by COs» in the atmosphere 
is important in meteorology. Unfortunately our 
measurements on CO, are less reliable than those 
on other gases because of the difficulties involved 
in removing atmospheric CO, from the optical 
path of our spectrograph. 

The CO: samples were admitted directly to the 
absorption cell from a Matheson lecture bottle. 


INTEGRATED APPARENT 
ABSORPTION iN 


CYCLES 
2350 CM BAND 
200 
150 
670 CM BAND 
FROM DATA OF 
100 + MARTIN AND BARKER 
670 CM BAND 
EQUIV PATH LENGTH IN CM AT N.T.P 
05 10 FOR 2350 cM BAND 


12 FOR 670 RAND 


Fic. 1. Experimental data on absorption by carbon dioxide 
bands at a total pressure of one atmosphere. 


3 These are summarized by O. Fuchs, Zeits. f. Physik 
46, 519 (1927). 
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INTEGRATED APPARENT 
ABSORPTION 
CYCLES «10 


200 F cm” BAND 


150 670 CM BAND 
100 
50 
Equiv, PATH LENGTH IN CM AT NTR 
75 70 FOR 2380 cm 
i iL. iL 
A 2 3 4 5 -6 FOR 670 CM BAND 


F 1G. 2. Experimental data on absorption by carbon dioxide 
bands at a total pressure of four atmospheres. 


The broadening gas was nitrogen obtained from 
an Air Reduction Company high pressure 
cylinder. One set of measurements was made at a 
total pressure of one atmosphere; a second set at 
about four atmospheres total pressure using a 


pressure cell. One atmosphere was not enough 


to give complete pressure broadening, but was 
probably sufficient to give a fairly close approach 
to the limiting absorption curve which would 
have been obtained with complete broadening. 

It is obviously necessary to eliminate the ab- 
sorption by atmospheric CO, in the optical path 
of the spectrograph in order to measure the ab- 
sorption of the CO, sample in the absorption cell 
even when a blank comparison cell is used. The 
presence of atmospheric COz2 introduces an inad- 
missable frequency dependence of the incident 
intensity.! The atmospheric COz would, in fact, 
absorb quite completely all frequencies corre- 
sponding to the lines in the CO». bands, so that 
addition of CO: in the absorption cell would 
have very little effect upon the transmitted light. 
The presence of atmospheric CO» would, there- 
fore, make the absorption by the CO, in the 
sample cell much smaller than it should be. 

In order to remove the CO: from the optical 
path, the spectrograph? was enclosed in a housing 
built of Masonite and plywood.*** The air inside 
was circulated through Drierite and Caroxite by 
means of a small blower. At the same time 
nitrogen was admitted to the housing at a rate 
of about 10 cubic feet per hour so as to maintain 
an outward flow through the many small leaks. 


*** Dr. Oliver A. Schaeffer assisted in the setting up of 
this air-conditioning system. 
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TABLE I. Intensities of carbon dioxide bands 
(units are cycles per cm at N.T.P.). 


670 2350 

From data of Martin and Barker® 660 X 1010 —— 
From our data with total pressure: 

1 atmosphere 560 X 10” 7200 X 1018 

4 atmospheres 470 X10 9600 x 10 

4 atmospheres 8900 X 
Average 560 X 1010 8600 « 10” 
From dispersion data of Fuchs¢ 538 X10” 7680 X10” 


* See reference 4. 

b . second set of data on the 2350 cm~! band was taken giving this 
result. 

© See reference 3. 


This procedure, while inconvenient, reduced the 
CO: concentration so that no absorption by at- 
mospheric CO: could be observed. In the first set 
of measurements, however, it was not possible to 
eliminate the atmospheric absorption com- 
pletely, and it was necessary to introduce a cor- 
rection for what remained. 

The experimental data are given in Fig. 1 and 
Fig. 2. In Fig. 1 we also show a curve based on 
data published by Martin and Barker‘ for the 
670 cm™ band. Their resolving power was sub- 
stantially greater than ours. Since the band at 
670 cm™ has a strong Q-branch, ‘envelope 
effect” is to be expected,! and it is likely that the 
difference between the two curves is due to the 
difference in resolving power. 

The data given in Fig. 1 and Fig. 2 are values 
of the “integrated apparent absorption” defined 
by Eq. (4), where ®=“‘integrated apparent ab- 
sorption,”’ 


@= f In(To/T)dy, (4) 


4 band 


To=light transmitted by blank comparison cell, 
T=light transmitted by sample cell. The ‘‘in- 
tegrated apparent absorption coefficient,” B, 
defined by Eq. (5) 


B=@8/pL (5) 


was calculated, as before,? from points on the 
smooth curves. 

The resulting intensities are given in Table I. 
The different values are not in particularly good 
agreement. The experimental errors may, how- 
ever, be somewhat larger for COz than the 10 
percent previously estimated,? so that the dis- 


4P. E. Martin and E. F. Barker, Phys. Rev. 41, 291 
(1932). 
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crepancies are no larger than should be expected 
on the basis of experimental errors. 

In the case of the 670 cm= band, difference 
bands are listed as occurring at 618, 647, 668, 
720, and 740 cm—.® We have not, however, tried 
to correct the observed intensity for contribu- 
tions from these weak overlapping bands. 

The comparison between our values of the 
intensity and those calculated from Fuchs’ 
expression for the index of refraction is made in 
Table I. The agreement is satisfactory, con- 
sidering the larger experimental errors for this 
substance. 

If our intensities are used to calculate the 
atomic polarization, we obtain 0.68 cc. The value 
given by Watson and Ramaswamy’ is 0.81 cc. 
Again the agreement is about as good as we 
expect. 

The properties of the CO bond moment may 
readily be calculated from the observed inten- 
sities. Since the bending and stretching vibra- 
tions have different symmetry, it is only neces- 
sary to calculate the effective masses for the two 
modes. The final values obtained for the bond 
dipole moment are: du’/dr 
= +6.010-'. These are similar to the results 
obtained for N,O.? 

This result for nu’ may be compared with the 
conventional values of 2.310-'8 assigned to 
C=O, and 0.7X10-' assigned to C—O.’? Ac- 
cording to Pauling,’ resonance in CO, involves 
three structures in about equal proportions: 


A B C 
0=C=0 6=C-0 0-C=6 


Since the formal charges cancel out, the observed 
value of u’ would be the average of that for 
single, double, and triple bonds. Our low value 
for wu’ implies a small dipole moment of the 
triple bond. It is evident from the three resonant 
structures for CO, that stretching a CO bond 
will tend to make it a single bond (which has the 
greatest length), negative at the oxygen end, so 


5 See, for example, G. Herzberg, Infrared and Raman 
Spectra of Polyatomic Molecules (D. Van Nostrand Com- 
pany, Inc., New York, 1945), p. 274. 

6 H. E. Watson and K. L. Ramaswamy, Proc. Roy. Soc. 
(London) A156, 144 (1936). 

7S. Glasstone, Textbook of Physical Chemistry (D. Van 
Nostrand Company, Inc., New York, 1940), p. 542. 

8L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, Ithaca, 1945), second edition, p. 196. 
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that du’/dr should be large, as it is, and negative. 
Since oxygen in more electronegative than 
carbon, we should select negative signs for both 
and dp’ /dr. 


METHANE 


In any study of the properties of hydrocarbons 
methane is a logical starting place. The inten- 
sities of the methane bands were not, however, 
the first to be measured, since it was anticipated 
that it would be very difficult to secure adequate 
pressure broadening for methane. The individual 
rotational lines would be widely spaced because 
of the molecule’s small moment of inertia, and a 
pressure of one atmosphere would not be ex- 
pected to provide sufficient pressure broadening. 
With the pressure cell, however, we were able to 
employ pressures up to five atmospheres. By 
using ammonia as a broadening gas considerably 
more effective pressure broadening was achieved 
than in the case of non-polar gases. In this way 
it was possible to make satisfactory measure- 
ments of intensities with methane. 

The methane was used directly from a Mathe- 
son lecture bottle (C.P. grade), without further 
purification. Ammonia used for pressure broad- 
ening was obtained from the Ohio Chemical and 
Mfg. Company. The experimental data are 
plotted in Fig. 3. Values of the ‘integrated ap- 
parent absorption coefficient’? calculated from 
these curves lead to the intensities given in 
Table II. 

In Table II a comparison is again made with 
data based upon measurements of infra-red dis- 
persion. The intensities obtained by the two 
methods are in good agreement, the differences 
being well within our experimental error. 

As a final check upon these intensities for 
methane, we may calculate the atomic polariza- 
tion from them. The result is 0.091 cc, whereas 
Watson and Ramaswamy’ obtained 0.08 cc. The 
agreement is good. 

A normal coordinate treatment of methane 
was made by Rollefson and Havens? to determine 
values of uw’ and dy’/dr for the CH bond. Our 
values are only slightly different, and were ob- 
tained by adjusting the values of Rollefson and 
Havens in proportion to the square root of the 


9 R. Rollefson and R. Havens, Phys. Rev. 57, 710 (1940). 
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INTEGRATED APPARENT 
300 | ABSORPTION IN 
cycles 
250 3020 CM BAND 
200 + 
150 |} 1306 BAND 
100 F 
50 
A 
EQUIV. PATH LENGTH IN CM AT NTR 


Fic. 3. Experimental data on absorption by methane bands. 


intensity. The results are up’ =0.31 K10-'8, du’ /dr 
= +0.55 X10-". 

These results are in fair agreement with the 
two values previously reported for the CH bond 
in ethylene, namely, 0.37 X 10-8 and 0.52 10-"8 
for uw’, about +0.610-' for du’/dr.2 As before, 
the sign of dy’/dr is undetermined. We have 
followed the usual convention in assuming the 
bond moment, y’, to be positive at the hydrogen 
end. T heoretical calculations by Coulson!® led, 
however, to a negative value for yu’. Our results 
do not give any information concerning the sign 
of p’. 


ETHANE 


* Measurements were made on ethane following 
those on methane, since it was desired to compare 
the behavior of the CH bond in the two related 
molecules. Ethane was used from a Matheson 
lecture bottle without subsequent purification. 
One atmosphere of air was adequate as broaden- 
ing gas. No attempt was made to reduce atmos- 
pheric absorption, which was appreciable only 
for the bands near 1400 cm~. 

In these experiments it was not possible to 


TABLE IT. Intensities of methane bands 
(units are cycles per cm at N.T.P.). 


Intensity Intensity from 
from these dispersion 

Band measurements measurements* 
1306 450 x 101° 444 x 10" 
3020 900 101° 970 10" 


* See reference 9. 


10 C, A. Coulson, Trans. Faraday Soc. 38, 433 (1937). 
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tie FF intense in ethylene but rather weak in ethane. 
z . . 
ca a Thus the abnormally great intensity of the out- 
of-plane bending mode of ethylene" has disap- 
gs peared in ethane. The stretching modes, how- 
one | Be ever, appear to be somewhat more intense for 


2950 cM” BAND 


3 


° 


1950 CM BAND 


200} So} $20 CM BAND 


EQUIV, PATH LENGTH IN CM AT NTR 
5 LS 


Fic. 4. Experimental data on absorption by ethane bands. 


resolve the band at 2954 cm from that at 
2994 cm=, nor the band at 1379 cm from that 
at 1486 cm™. Each of these overlapping pairs of 
bands is, therefore, considered as a single band 
in tabulations of experimental data. The fre- 
quencies are listed as 2950 cm~ and 1450 cm“, 
respectively. The data are presented graphically 
in Fig. 4. Values of the “integrated apparent 
absorption coefficient” calculated from these 
curves lead to the intensities given in Table III. 

Unfortunately there are no other data on 
ethane intensities available which could be used 
to confirm or deny these results. A very rough 
check is possible in terms of atomic polarization. 
Our intensities lead to a calculated atomic 
polarization of 0.14 cc, whereas the value given 
by Watson and Ramaswamy’ is 0.09 cc. This is 
fairly good agreement, since small atomic polari- 
zations could not be measured very accurately 
by Watson and Ramaswamy’s method. 

Values of the intensities of ethylene are also 
given in Table III. The most striking difference 
is in the low frequency bending mode which is 


TABLE III. Intensities of ethane bands 
(units in cycles per cm at N.T.P.). 


Band Intensity 
Ethane: 
820 73X10! 
1450 230 X 
2950 cm~! 2200 x 
Ethylene :* 
950 1540 10" 
1444 cm“ 210 10" 
3050 cm= 840 x 101° 


* From reference 2. 


ethane than ethylene. 

A more detailed interpretation of ethane in- 
tensities in terms of properties of the CH bond 
can be made if the intensities of the 1450 cm 
and 2950 are broken down into 
their individual components. Unfortunately, no 
good basis for estimating the intensities of the 
individual overlapping bands is available. The 
best published data appear to be those of Levin 
and Meyer” who, however, failed to resolve 
either pair completely. Qualitative estimates of 
the intensities are given by Herzberg." 

In the case of the 1379 cm™ and 1486 cm 
bands his estimates are ‘‘weak’’ and “strong,” 
respectively. Levin and Meyer’s data indicate 
fairly clearly that the 1379 cm~ band is a good 
deal the weaker, and our absorption curves show 
a weak maximum at about 1375 cm. As a 
rough estimate, therefore, we have assigned 
30X10! to the 1379 band and 20010!" 
to the 1486 cm™ band. 

In the case of the 2954 cm™ and 2994 cm 
bands Herzberg’s designations are “‘strong’’ and 
“medium” respectively. The experimental data 
of Levin and Meyer are consistent with this 
estimate, but could apparently equally well 
arise from a weak band at 2954 and a strong one 
at 2994. In our data the two are completely 
unresolved. Qualitative consideration of the 
geometry of the molecule suggests, however, that 
the band at 2994 cm is the stronger, since it is 
a perpendicular band.'* On the basis of the 


11 Data on ethylene are discussed in detail in reference 2. 
12 A, Levin and C. F. Meyer, J. Opt. Soc. Am. 16, 137 
(1928). Measurements at higher resolution have been 
reported (L. G. Smith and W. M. Woodward, Phys. Rev. 
61, 386 (1942)) but no experimental curves given. 

13 See reference 5, p. 344. 

14 The reasoning is as follows: In the ethane molecule 
the CH bonds are almost perpendicular to the CC bond 
(the C—C—H angle being about 109°). Therefore a CH 
bond stretch will give rise to a contribution to uw whose 
component parallel to the CC bond is small and perpen- 
dicular to it is large. In a parallel type band the perpen- 
dicular components cancel out, while the parallel compo- 
nents add up to a small net dipole moment. In a 
perpendicular-type band the parallel components cancel 
out, while the perpendicular components add up to a large 
net dipole moment. Therefore the perpendicular type band 
should be the stronger one. 
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TABLE IV. Apportionment of ethane intensities. 


Frequency Intensity in cycles 
Band in cm~! per cm at N.T.P. du/3Q 
V5 2954 200 x 10%” + 45.9 
VE 1379 30 10” + 17.8 
v7 2994 {2000 x 10” +103 
vs 1486 20010" + 32.5 
Vy 820 { 73x10" + 19.7 


normal coordinate treatment, which will now be 
discussed, it was determined that the most con- 
sistent results were obtained by assigning an 
intensity of about 200X10'° to the 2954 cm 
band and 2000X10!° to the 2994 cm- band. It 
should be emphasized, however, that these inten- 
sities are tentative estimates; only their sum has 
been accessible to direct measurement. 

A simple normal coordinate treatment was 
carried out for ethane in order to determine the 
properties of the CH bond moment. It was 
assumed, following Howard,'* that the potential 
function involves only two force constants and 
may be written as in Eq. (6). 


2V=K, >, B?+> |. (6) 


jk 


Here r;=change in CH; bond distance, B; 
=change in C—C—H,; angle, and a;,=change 
in H;—C—H, angle. In the latter case the two 
hydrogens must obviously be attached to the 
same carbon. K,;=4.79X105 dyne/cm, K2=0.46 
X10° dyne/cm. The method and terminology 
used follow that of Wilson.'® 

The configuration of the ethane molecule is 
uncertain to the extent that it may be either 
“staggered” or ‘‘eclipsed.’’ Both lead, however, 
to the same relations between the observed inten- 
sities and properties of bond moments.!” Our cal- 


16 J. B. Howard, J. Chem. Phys. 5, 442 (1937). 

16 FE. Bright Wilson, Jr., J. Chem. Phys. 9, 76 (1941). 

17 This may be shown in the following way: We can 
change from staggered to eclipsed forms by rotating one 
CH; group by 180° about the CC axis. The symmetry 
coordinates for vibrations parallel to this axis are un- 
changed, those for vibrations perpendicular to the axis 
must have coefficients of the internal coordinates belonging 
to the rotated CH; group changed in sign. For example, 
the symmetry coordinate @za is given by 1/v3[r—r;s 
—4(r2+r3—rs—re) | for the staggered case as contrasted 
with 1/v3[71+74—43(ro+7r3+75+7¢) ] for the eclipsed case. 
In order to set up the secular equation, following the 
method of reference 16, we calculate the vectors, S defined 
by Eq. (11) of reference 16 as S;%=2Z,% UiSiz. In the 
case of Ra the s vectors are unit vectors along the CH 
bonds. The Ux, are coefficients which give Rr in terms 
of the r’s. When one changes from eclipsed to staggered 
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TABLE V. Characteristics of CH bond moment in ethane 
(c.g.s. units). 


Type of 
vibration Op’ /dr 
ve-bending 0.104 X 10718 or 0.096 or 
0.094 10718 0.086 X 10-1 
vs-bending 0.35 10718 or 0.32 o1 - 
0.33 0.30 
vy-bending 0.30 or 0.28 or 
0.26 0.24 
vp-stretching +0.69 107" or 
—0.74 
v7-Stretching -+0.81 X10-” or 
—0.74 X10-" 


culations were made on the basis of the eclipsed 
configuration, but the final results are equally 
valid for the staggered case. 

The symmetry coordinates for the infra-red 
active vibrations were taken to be: 


Parallel-type modes 
Re = —Bs —Bo 


— (a23+ —a56 — — 045) |. 


Perpendicular-type modes 
Rre = 
= (374) +Bs— 
(Rso= 
Roa = (374) |, 


= 3 (ais —as5). 


In terms of these symmetry coordinates the 
normal coordinates are given by the trans- 
formation L~! whose elements are tabulated 
below : 


(L™)55= —1.276XK10-", —0.0177 
(L™)77= 1.225 X (L)73 = —0.0040 10-, 
(L~)79=0.0237 X 10-™, 
(L~)s7= —0.0355X10-", (L~)ss=0.198X 
(L~*)s9= 0.808 X 10-, 
—1.313X 
(L)99=0.322 X 10-. 


(L~)o7 =0.0914 X 10-”, 


configuration some of the s vectors change direction, but 
the corresponding Ui, coefficients change sign so that there 
is no net change in the S vectors or in the secular equation 
calculated from them. Either form for the molecule will 
then lead to the same set of L~' coefficients (in terms of 
symmetry coordinates), and the same values of the 
(0u/dRx)’s determined from the observed intensities. 
The interpretation of the (0u/d®,)'s in terms of the 
(0u/dR)’s, which are properties of individual bonds, 
again involves the geometry of the molecule and the Ux, 
coefficients, but again in such a way that the changes in 
these two factors compensate in going from eclipsed to 
staggered form. Therefore the same results will be obtained 
for properties of bond moments whichever form of molecule 
is assumed. 
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The procedure for calculating the properties 
of bond dipole moments is the same as that 
previously described.? The perpendicular modes, 
v7, vg, and vy are each twofold degenerate, so 
that only half the observed intensity should be 
ascribed to each member of these degenerate 
pairs. The resulting values of intensities and of 
(du/dQ)’s are given in Table IV, (du/0Q)’s 
being calculated by Eq. (3). 

The values of the (du/0Q)’s and (L~")’s aré 
substituted in Eq. (7) in order to obtain values 
of the (du/dR)’s. From these 


dy = x (7) 


latter and the geometry of the molecule the 
dipole moment, yw’, and its rate of change with 
internuclear distance, du’/dr, may be calculated. 
If we exclude negative values for the bond dipole 
moment as unreasonable, the final values are 
those given in Table V. (Here s is the equilibrium 
CH distance.) If du’/dr is positive, the choice of 
signs for the (du/0Q)’s must be that leading to 
the top member of each pair of values in Table V. 
If negative, the bottom member is to be taken. 

These results may be compared with those 
obtained for methane. The outstanding difference 
between the two molecules is the very low value of 
un’ obtained for the parallel-type band v¢ of ethane. 
The intensity of this band is only about #g what 
would be required to give the expected value of yu’. 
The reason for this is not known and merits 
further investigation.!® The values of yu’ obtained 
from vg and vg agree well with that from methane, 
however, all being close to. 0.3X10~'*%. In the 
case of du’ /dr the agreement is less satisfactory ; 
values from ethane are about +0.75X10~-", 
whereas +0.55X10-! had been obtained for 
methane. 

These results cast some doubt upon the whole 
concept of additivity of bond dipole moments. 


18 Such a large discrepancy as this could not be accounted 
for as an error in dividing the intensity of the ‘1450 cm™ 
band” between ‘vg and vs. To give vg ten times as much 
intensity we would have to assign to it all the intensity of 
the “1450 cm™ band” and more, leaving nothing at all 
for vs. 
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It would seem that the apparent values of yu’ and 
Ou’/dr change so much from one molecule to 
another and one type of vibration to another 
that strong interactions with neighboring bonds 
must frequently occur. One might think, how- 
ever, that our normal coordinate treatment of 
ethane has been too simple. The potential func- 
tion with only two force constants does not give 
very good solutions for the secular equation, and 
may introduce very considerable approximations 
in the normal coordinates. Nevertheless the fact 
remains that the stretching modes are, relatively, 
too intense in ethane, and the bending modes too 
weak. The conclusions concerning p’ and dy’/dr 
are actually not at all sensitive to changes in the 
normal coordinate treatment. 


CONCLUSIONS 


Two fairly definite conclusions may be drawn 
from the results summarized in this and the pre- 
ceding paper: 

(1) Our intensities, obtained by direct meas- 
urement of absorption employing an extrapola- 
tion to zero absorption, are in satisfactory agree- 
ment with results obtained from measurement 
of infra-red dispersion in all cases where a com- 
parison is possible. Our method may, therefore, 
be considered reliable and applied with con- 
fidence to other molecules. 

(2) It is not yet possible to predict intensities 
of vibrational bands on the basis of empirically- 
assigned dipole moments for the chemical bonds 
making up a molecule. If this could be done, it 
would greatly extend the usefulness of infra-red 
spectroscopy in the analysis of molecular struc- 
ture, since the infra-red intensities could be used 
to determine the numbers of various types of 
structure present. Further investigations are 
necessary to measure and explain the behavior 
of a variety of bond dipole moments under a 
variety of conditions. 

In conclusion the author takes pleasure in 
acknowledging the continued help of Professor 
E. Bright Wilson, Jr., under whose guidance 
these investigations were carried out. 
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ail A Note on the Relation between Entropy and Enthalpy of Solution 
to O. K. Rice* 
1er Clinton Laboratories, Oak Ridge, Tennessee 
ids (Received June 2, 1947) 
- 
of The relation between entropy and enthalpy of solution for a series of non-polar solutes in a 


given non-polar solvent is discussed. It is considered that solution of a gaseous solute, without 
changing its concentration on going from gas phase to liquid solution phase, does not change 


ive its own entropy, all change of entropy being referred to the solvent. The entropy of the solvent 
ind changes because of surface effects around the solute molecules and a kind of long-range order 
ons introduced by the solute molecules. Two extreme cases are considered, (1) the case of an ideal 
act solution, and (2) the case of a solute of hard attractionless spheres. The difference in entropy 
rl of solution between these extreme cases can be estimated. It can also be estimated by extra- 
ny polation from the experimental data on the entropy and enthalpy of solution, and these two 
too estimates agree in order of magnitude. The fact that the relation between entropy and enthalpy 
/dr of solution is linear is also shown to be a reasonable expectation, and effect of changing solvent 
the as well as solute is considered. The groundwork is thus laid for a qualitative understanding of 
this relation between entropy and enthalpy of solution. 
iia T is well known that there is a relation between _ being the case, it appears that some consideration 
cam the entropies and enthalpies of solution of a of this problem from another point of view is 
series of non-polar gases or vapors in a given _ still desirable. 
a non-polar liquid solvent, provided the solute In the present note we shall consider the 
“a molecule is not larger than that of the solvent. properties of two extreme, hypothetical types of 
sai This was apparently first remarked by Bell,| who solution, namely (1) an ideal solution, and (2) a 
int noted that some relation between the energy and _ solution in which the solute molecules consist of 
-_ enthalpy of solution might be expected but gave hard spheres, which occupy space but do not 
re, no explanation of the form of the relation nor of exert attractive forces on the solvent molecules. 
on. the order of magnitude of the quantities in- Real solutions may usually be considered to 
volved; it was further discussed by Barclay and stand somewhere between these extreme types, 
on Butler.2 A rather detailed theoretical discussion and it is possible to say something about the 
lly- has been given by Frank and Evans.’ In their transition from one type to the other for any 
seth treatment certain relations used in the discussion given solvent. 
2, It of pure liquids? bdaesind introduced into the theory Notation.—Since we shall be interested in entropies of 
_red of solutions only slightly, if at all, modified.5 It gas and liquid phases for solvent and solute and changes 
na is certainly of interest that this gives an equation of entropy under various conditions, as well as a number 
iil which seems to fit pretty well the energy- of different preps and free volumes, we require od system 
| sitadins Sey citations acted chews: of notation with a variety of subscripts and suffixes, and 
s of enthalpy - Shh .,. capital and sinall letters. In the table which follows, a 
are however, the physical significance and validity subscript i is used to denote component i. In the text i 
vior of this procedure are not entirely clear. This takes the value 1 for solvent and 2 for solute. To any of 
isa —entenaei the symbols we may affix an asterisk, indicating that it 
; *On leave, 1946-47, from the University of North — refers to the case of an ideal solution, or we may add a 
Carolina, Chapel Hill, North Carolina. This document is subscript 0, indicating that it refers to the case in which 
> in based on work performed under Contract No. W-35-058- the solute consists of hard-sphere molecules. 
eng-71 for the Atomic Energy Project at the Clinton 
ssor Laboratories. n,=total number of moles of component i. 
ance 1 R. P. Bell, Trans. Faraday Soc. 33, 496 (1937). N;=mole fraction of component 7 in liquid phase. 


2]. M. Barclay and J. A. V. Butler, Trans. Faraday 
Soc. 34, 1445 (1938). 

3H. S. Frank and M. W. Evans, J. Chem. Phys. 13, 
507 (1945). 

4H. S. Frank, J. Chem. Phys. 13, 493 (1945). 

5 For example, the use of Eq. I1(10), together with the 
subsequent use of 8; and £2 as partial molal quantities on 
p. 508 of Frank and Evans; and the use of Eq. II(21) on 
p. 513. 


V=volume of liquid solution containing one mole of 
solute. 

Vi°=molal volume of pure liquid component i. 

v;=molal volume of gaseous component 7. 

Vs; =molal free volume of component 7 in liquid solution. 

Vy,i° =molal free volume of component i in pure liquid 
phase (see Eq. (5)). 
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S=total entropy of liquid solution. 

8;°=molal entropy of pure liquid component 7. 

8;=partial molal entropy of component 7 in liquid solu- 
tion. 

si(¥;)=molal entropy of gas at volume indicated. 

AS;(v;) =8; —s;(v;) =entropy of solution of 1 mole of pure 
vapor component 7 at molal volume v; in a large 
amount of solution of any given concentration 
(partial entropy of solution). 

—si(v;) =molal entropy of condensation of 
pure component 7 from molal volume 7; in vapor 
phase. 

ASm=entropy of mixing. 

AH; =enthalpy of solution of 1 mole of pure vapor com- 
ponent 7 in a large amount of any given solution 
(partial enthalpy of solution). 

AH;°=enthalpy of condensation of pure component 7. 


IDEAL SOLUTIONS 


If two liquids form an ideal solution, the 
entropy of mixing is given by the classical 
expression 


= InNi—72.R InNa, (1) 
and the total entropy of the solution is given by 
S* InNi—72R InNo. (2) 


We then obtain for the partial molal entropy of 
the solute 


= (0S*/dn2) my =S.°—R InNo. (3) 


Let us now consider the entropy of solution of 
1 mole of solute from the vapor phase into a 
very large quantity of solution of any given 
concentration, with v;= V. This means that the 
solute enters the solution without change in its 
volume concentration. We have for this process 


AS2*( V) =§,* —So( V) 
= S.° — So( V) —R InNo. (4) 


We can write an equation defining V; 2° as 


follows: 
So° — Se(v2) = R In(Vy, 2°/V2). (5) 


The free volume thus defined includes all effects 
of communal entropy. 

An exactly ideal solution is one in which the 
solute and solvent are ‘exactly alike” although 
distinguishable ; that is, their molecules have the 
same size and force fields.* In this case = V2? 


6 There is, of course, no requirement from the thermo- 
dynamic definition of an ideal solution (giving the chemical 
potentials as functions of the mole fractions) that the 
molecules of solvent and solute be ‘‘exactly alike,” and 
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and V;,1°=Vy,2°. Since V is the volume of solu- 
tion containing 1 mole of solute, we have 
V=V.2°/Ne. Setting ve= V, introducing the ex- 
pression for V into Eq. (5), and substituting in 
kq. (4), we obtain 


AS)*(V) = R In(Vy,2°/V2°) 
=R In(Vy, 1°/Vi1°). (6) 


The ratio of molal free volume to molal 
volume simply represents the effect of neighbor- 
ing molecules on the space available for the 
motion of any molecule in the field of those 
neighbors. Thus the effect of solution without 
change of concentration is simply the effect of 
neighboring molecules on the free space available. 


NON-IDEAL SOLUTIONS 


Equation (6) holds, of course, only for ideal 
solutions. We may, however, define a new 
quantity V,,;’ by the general equation 


AS.(V)=R In(V;, 1’/V1°). (7) 


V;,1’, then, may be said to give that volume left 
free by 1 mole of solvent for a particular solute. 
It seems possible, as Frank and Evans have 
noted, for V;,:’ to become as large as, or even 
greater than, Vi°. At first sight this appears 
strange, for it seems then that the solute is free 
to move around in a greater volume than that 
which contains it, even though this volume itself 
is well filled with solvent molecules. Frank and 
Evans pointed out that this could only be 
explained as an effect of the solute on the solvent. 
An explanation of the nature of this effect, and 
an estimate of its order of magnitude, will be 
the principal aim of the present note. 

We now consider the process of solution from 
a different point of view. We suppose that we 
have the solute in the volume which it is going 
to occupy, and pour the solvent in on it. Since 
this is a liquid system, dilute in the solute, the 


there are many solutions, that are very close to ideal, for 
which this is certainly not the case. However, for a solution 
to be quantitatively ideal, it seems highly probable that 
this condition is necessary. It certainly seems to be 
necessary as far as the size of the molecules is concerned, 
though this is not very critical. [(See R. H. Fowler and 
E. A. Guggenheim, Statistical Thermodynamics (Cambridge 
University Press, 1939), pp. 353 ff.) ] The effect of the 
force fields on the ideality of a solution as expressed by 
Eq. (1) is essentially the subject of the present note. We 
shall, hereafter, refer to the “exactly ideal” solution 
described here simply as an ideal solution. 


| 
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TABLE I. Enthalpies and entropies of solution with acetone as solvent. 


Hard-sphere solute Ideal solution 


Partial enthalpy of solution AH2,o= 


Partial entropy of solution 
Total entropy of solution 


—AH,° = 7600 
AS2,o( V) =8 AS:*(V) = —11.5 
V)+R=10 


=AH,° = —7600 


Enthalpies in cal. mole, entropies in cal. 


solute molecules may occupy any preassigned 
positions, regardless of whether the solvent is 
present or not. We may thus say that the parti- 
tion function of the solute is unaltered by the 
presence of the solvent, and we may, somewhat 
artifically perhaps, refer all entropy effects to the 
solvent.’ In considering the entropy of the sol- 
vent, we may assume that all solute molecules 
are held in fixed positions, since in the partition 
function for the solute all possible positions are 
included. From this point of view Eq. (7) gives 
the effect of the solute on the entropy of the 
solvent; however, a slight correction is required. 
For it is evident that the process we have just 
considered gives the total entropy of solution of 
1 mole of solute rather than the partial entropy 
of solution. The total entropy change will be 
S—s2(V)—8:° instead of AS2(V)=S.—s2(V). 
Since (since m2=1), we 
see that the entropy change for the process 
considered is equal to AS2(V)+-8:—,8,°. For 
a dilute solution this may be shown’ to be equal 
to AS2(V)+R. 

For a solution that is ideal as well as dilute, 
the total entropy of solution AS,*(V)+R of one 
mole of solute molecules will be negative; this 
follows from Eq. (6) because V;,1° is always very 
considerably smaller than V,°. But in an ideal 
solution the solvent molecules in the immediate 
neighborhood of a solute molecule are in the 
same environment as any other solvent mole- 
cules, since the force fields are the same. The 
lowering of entropy of the solvent implied by 
Eq. (6) cannot, therefore, be referred to any 
change in the range of motion of the individual 
solvent molecules in the neighborhood of solute 


7 Somewhat similar considerations have been carried out 
by R. H. Fowler and E. A. Guggenheim, pp. 372 ff of 
reference 6, but they did not attempt to consider the 
solvent effects fully in setting up the chemical potentials. 
See also Barclay and Butler, reference 2, p. 1454. 

8 We have §;=8:°—R InN, since Raoult’s law holds for 
a dilute solution. Ny=;/(m1+1) with n.=1. With 
this becomes and InN; ~—m,~!, whence the 
relation follows immediately. 


molecules. It must, on the contrary, be attributed 
to the introduction of a certain degree of long- 
range order, produced by having the solute 
molecules held in fixed position. This restriction 
in position is transmitted through the neighbor- 
ing molecules to the solvent, even though the 
range of motion of these molecules about their 
equilibrium positions is not altered. It is rather 
the equilibrium positions themselves which are 
affected. 

To understand better the situation in non-ideal 
solutions, let us consider the case in which the 
solute is assumed to be composed of hard-sphere 
molecules, which exert no force on the solvent 
molecules. (This means that the enthalpy of 
solution, AH>2» is actually positive because of 
the energy necessary to produce the hole in the 
solvent into which the solute molecule is going 
to go. Roughly, assuming that the hole is the 
same size as a solvent molecule, we may say 
that A/T» is equal to —AH,°.)° If the solute is 
like a hard sphere, we may expect that a fixed 
arrangement of solute molecules will be much 
less effective, if effective at all, in inducing long 
range order in the liquid. Furthermore, since the 
energy will not be lowered by proximity of 
solvent to solute molecules, the solvent mole- 
cules around a solute molecule will be reasonably 
free to arrange themselves in such a way as to 
allow a maximum of freedom of motion. There 
will thus be a gain in entropy, similar to the gain 
in entropy when the free surface of a liquid is 
increased. This gain in entropy is to be equated 


to o(V)+R. 


APPLICATION TO THE EXPERIMENTAL DATA 


Let us now apply these ideas to dilute solutions 
in a typical solvent, acetone. We use the data 
collected by Frank and by Frank and Evans,’° 


®For a calculation of the energy necessary to form a 
hole in a liquid, see R. H. Fowler, Statistical Mechanics 
(Cambridge University Press, 1936), second edition, p. 844. 

10H. S. Frank, ref. 4; see especially p. 499; H. S. Frank 
and M. W. Evans, ref. 3; see pp. 514-515. 


“a 
lal 
or- 
he 
se 
ut 
of | 
le. 
7) 
oft 
Ee. 
ve 
en 
urs 
ee 
at 
elf 
nd 
be 
at. 
nd 
be 
ym 
ve 
ng 
ce 
he 

for 
ion 
lat 
ed, 
nd 
lge _ 
the 
by 
ion. 


878 0. K. 


and present the relevant material in Table 1. 
The enthalpy values are obtained as indicated 
from the enthalpy of condensation of acetone. 
To get AS,*(V) we take the entropy of vaporiza- 
tion and the density, and using Eq. (5) (applied, 
however, to the solvent instead of the solute) 
calculate that at 25°C the value of V;,1°/V,° is 
0.0030. We then use Eq. (6). The table of 
Frank and Evans shows how AS.(V) and AH, 
vary from solute to solute. (Actually they list 
the energy and enthalpy for evaporation from 
solution to form a gas at 1 atmosphere at 25°C, 
so their values differ by a sign and a constant 
additive amount from ours.) By extrapolation 
we find that the value of AS» o(V) given in the 
table is that which corresponds to 


AH, +7600. 


Had there been no surface effect AS2o(V)+R 
would have been zero, assuming that there is no 
long-range order under these conditions. It is, 
therefore, natural to compare the 8+ R=10cal. 
mole deg.—! with the surface entropy to be ex- 
pected. The surface entropy of a liquid is closely 
related to the Eétvés constant, —d(V,° *y)/dT, 
where vy is the surface tension. If we divide this 
by N!, where N is Avogadro’s number, the 
expression may be interpreted as the surface 
entropy per molecule. (This follows because y is 
the surface free energy per unit surface, and 
(Vi°/N)! is very close to the surface occupied by 
a molecule at the surface.) The Eétvés constant 
for acetone is, from data in the Landolt-Bérn- 
stein Tables, 1.8 erg mole~? deg.—', which makes 
— N-4d(V,1° equal to 2.5 X erg deg.—. 
Comparing this to the Boltzmann constant 
k=1.37X10-"* erg deg.—', we see that the surface 
entropy is a little less than 2R or 4 cal. deg. 
per mole of molecules at the surface. 

The surface area about a solute molecule is of 
the order of four times the area occupied by a 
molecule at the surface. One would not expect 
the entropy connected with unit area to be the 
same as for a flat surface, and the 10 cal. mole 
for the surface entropy in the solution of a hard 
sphere gas is certainly of the correct order of 
magnitude. 

We may look at this from a slightly different 
point of view. There will probably be about 10 
nearest neighbors about a solute molecule. Each 
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neighbor, therefore, has about 1.0 cal. deg. 
mole of surface entropy, about one-fourth as 
much as a free surface molecule. In a recent 
paper," we estimated the surface entropy per 
molecule of a liquid by a rough statistical calcu- 
lation and found a value about half as great as 
that given by the E6tvés constant. In our calcu- 
lation we neglected the decrease in density at 
the surface layer, which undoubtedly results 
in an increase in the surface entropy. On the 
other hand, this decrease in density is prob- 
ably not appreciable in the neighborhood of 
a solute molecule. Further, a molecule at the 
surface of a small spherical hole inside a liquid 
is not nearly as free as a molecule at a flat 
surface. Therefore, the value of 1.0 cal. deg. 
mole for the solvent molecules about the hard- 
sphere solute molecule seems entirely reasonable. 

We may now be in a position to understand 
the linear relation between AS,(V) and 
simply as the start of a series expansion. It has, 
as we have noted, been pointed out by Barclay 
and Butler that there should be some sort of 
relation between AS2(V) and AH, for a given 
solvent, and this is also obvious from the dis- 
cussion given above. We expand about the values 
for a hard-sphere solute, writing 


AS2(V) =AS», V) +a1(AH2— 9) 
(8) 


The question then reduces to a decision as to 
whether the first term in this expansion is the 
dominant one. AS:(V) can be divided into two 
parts, as discussed above, the surface entropy 
and the entropy (negative in sign, of course) 
associated with long-range order, and each one 
of these can be expanded in a series like Eq. (8). 

Let us consider the surface entropy first. As 
we have noted, this part is contributed by all 
the nearest neighbors of a solute molecule, and 
in the range from ideal solution to solution of 
hard spheres it goes from 0 to 1.0 cal. per mole 
per deg. for each nearest neighbor, in the case of 
acetone. This corresponds to a change in the 
effective free volume of each nearest neighbor 
molecule by a factor 1.65. Over this great a range 
the free volume change might be expected to 
deviate somewhat from being a linear function 
of the force exerted by the solute molecule on 

10, K. Rice, J. Chem. Phys. 15, 314 (1947). 
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its neighbors and hence” on AHsg, and the loga- 
rithm of the free volume, which determines 
AS2(V), would also be expected to deviate to 
some extent from being a linear function of the 
free volume itself. However, these deviations 
would not be expected to be exceedingly great 
even with a 65 percent change in free volume, 
and the experimental data do not actually cover 
more than about two-thirds of the range between 
the ideal solution and the hard-sphere solute. 
And the change in surface entropy contributes 
only a small part of the variation of AS2(V) 
between these extreme types of solution. The 
greater part of this variation is to be referred to 
the introduction of long-range order. But this is 
an effect which is actually shared among many 
molecules of solvent, so that the change in free 
volume for any one molecule will be so small 
that one need have no surprise if the experimental 
range does not extend beyond that in which the 
first term of the series expansion of Eq. (8) 
suffices. This, of course, is not a rigorous expla- 
nation of the linearity between and AS,(V), 
but does make it seem plausible. Frank and 
Evans obtained a linear relation, but this was 
done by carrying an empirical linear relation 
between the AH’s and AS’s of pure liquids into 
the equation for the solutions. 

The slope of the AS2(V) vs. AH, line, which 
may be expressed as 


a, =[AS*(V) — AS», o( (9) 


is almost the same for a considerable number of 
solvents. This may be understood in terms of 
the variation of AS.*(V) with AH,* for a series 
of different solvents. This, of course, is the same 
as the variation of AS.°(V) with AH.°. 

Let us write 6AS,*(V) and 6AH,* for the 
differences between the respective indicated 
quantities for two solvents. Then since AH»,¢ 
= —Al],* we will have 


6(AH* 0) 26AH2*. 


12 This depends upon the existence of a fairly simple 
relation between the force exerted by neighbor molecules 
and the mutual energy, which means that the potential 
energy curves must in all cases be of similar shape. 


This means that if a; is to be the same for the 
two different solvents, we should, from Eq. (9), 
have 


V) = AS», o( V)) 2a,6AH.*. 


5AS2,o(V) may be expected to be close to zero, 
because AS», o(V) is contributed entirely by the 
surface entropy, and the E6étvés constant has 
roughly the same value for most common sol- 
vents. We might thus expect to find 6AS,*(V) 
~2a,6AH,*. Actually, empirically, it is found 
that 6A.S.*(V)=a,6AH.*, since the slope of the 
AS,°(V) vs. AH,° line for different solvents is 
approximately the same as that of the AS,(V) 
vs. AH; line for solutions with a common solvent. 
But for the usual solvents 6AS,°(V) is such a 
small fraction of AS2*(V)—AS2o(V) that the 
empirical relation between AS,°(V) and AH,° 
does not require that the AS,(V) vs. AHz2 lines 
for the different solvents be appreciably different. 
It thus appears that the similar values of a; for 
the different solvents may merely reflect the 
similarity of the solvents used. 

On the other hand, the fact that the AS,°(V) 
vs. AH,° line has about the same slope as most 
of the AS.(V) vs. AH» lines for the various 
solvents is a remarkable fact, noted by Barclay 
and Butler,? but still not fully explained. In 
going from one solvent to another, for which 
AH,° is greater, there is no change in surface 
entropy, since this is zero with either pure solvent 
or ideal solution ; the decrease in entropy arising 
from long-range order, however, is caused by 
two factors, (1) the increased force exerted by 
the particular atom which is condensed into the 
solution, and (2) the increased force which all 
neighbor atoms exert on each other. On the 
other hand, in going from one solute to another 
which has a greater AH2 in the same solvent 
there is a decrease in entropy on account of the 
surface effect, but there is no decrease in entropy 
on account of factor (2) of the preceding sentence. 
The surface effect in the case of changing solutes 
must, therefore, approximately balance factor 
(2) in the case of changing from one pure solvent 
to another. 
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Halogen-Halogen Distances in Halogen- 
Polymer Complexes 


R. E. RUNDLE 
Depariment of Chemistry, Iowa State College, Ames, Iowa 
October 31, 1947 


ECENTLY West has found that many polymers form 

complexes with iodine, iodine bromide, and bromine.! 
Stretched fibers of these complexes give x-ray spacings 
characteristic of one-dimensional lattices with periods of 
3.10, 2.90, and 2.70A, respectively, characteristic of the 
halogen but essentially independent of the polymer. West 
has interpreted the one-dimensional spacings as inter- 
atomic distances in linear, polyhalogen chains. 

The single, covalent bond distances in Iz, I-Br, and Bre 
are 2.67, 2.47, and 2.28A, respectively,? far shorter than 
the spacings observed by West. If the above spacings are 
halogen-halogen interatomic distances, their bond num- 
bers are only about 0.17 for all three cases.* 

While it is easy to visualize a polyhalogen chain with 
links of bond order } (see below), it is hard to understand 
the virtual elimination of halogen-halogen bonds in the 
complexes as West’s suggestion demands. One possible 
explanation involves bonding of the halogen atoms to the 
polymer, but since the polymer fiber spacings are unre- 
lated to that of the halogen separation,' this is hardly 
possible. Moreover, polymer-halogen bonds are unde- 
sirable chemically, and will not explain many of the effects 
observed in the amylose (starch)-iodine complex.‘ 

Since, at best, only three orders of the lattice spacing 
have been observed,! it is by no means certain that the 
halogen atoms in the one-dimensional array are equally 
spaced. Hence, it is possible that halogen molecules I., 
I-Br, and Bre exist in the array, but with somewhat 
shortened van der Waals separations between molecules. 

It has been suggested that a dipole field of a polymer 
acting on a highly polarizable, linear array of halogen 
molecules might account for the stability of the halogen- 
polymer complexes.‘ Recently Mr. R. S. Stein and the 
author have examined this proposal more quantitatively,® 
and find that very reasonable fields due to the polymer 
will lead to very materially reduced van der Waals dis- 
tances between iodine molecules in the iodine complex. 
A period of 6.2A (twice the observed period) is quite 
acceptable for an iodine-iodine bonded distance plus a 
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non-bonded van der Waals distance. The proposal also 
accounts for the very low fugacity of iodine over the 
complex and the other unusual physical properties of these 
complexes. 

Though polyhalogen chains with equally spaced halogens 
are allowed by our proposal, they result from resonance of 
the bond between bonded and non-bonded neighbors as 
the van der Waals and bonded distances approach each 
other under the influence of the polymer field. A bond 
number of 3 would result for each bond, and halogen- 
halogen distances of 2.85, 2.65 and 2.46A are then ex- 
pected for the three complexes above. These distances are 
so far below observed spacings it is unlikely that the 
limiting case is approached in any known complexes. 

Our calculations allow variation in the halogen-halogen 
distance with dipole field of the polymer. West has noted 
variations in spacing of 0.15A. This variation is difficult to 
understand in terms of the polyhalogen suggestion. 

In summary, West’s observed lattice constants appear 
too long for polyhalogen chains, but are in excellent agree- 
ment with expectations based on the dipole interaction 
suggestion. Consequently, they tend to substantiate the 
latter proposal rather than otherwise. 

1C. D. West, J. Chem. Phys. 15, 689 (1947). 

2L. Pauling, Nature of the Chemical Bond (Cornell University Press, 
Ithaca, New York, 1945) second edition, p. 164. 

3 Calculations made in accordance with the rule R(x) =R(l) —0.3 
logn, as given by L. Pauling, J. Am. Chem. Soc. 69, 542 (1947). 

4 Rundle, Foster, and Baldwin, J. Am. Chem. Soc. 66, 2116 (1944). 


5 R. S. Stein and R. E. Rundle, ‘The nature of the interaction be- 
tween starch and iodine’ (to be published in J. Chem. Phys.). 


Infra-Red Dichroism Involving Hydrogen Bonds 


LEONARD GLATT AND JOSEPH W. ELLIs 
Physics Department, University of California, Los Angeles, California 
October 28, 1947 


T is well known that the vibration bands of O—H units 
which take part in hydrogen bonds are broadened and 
shifted toward longer wave-lengths. Most spectroscopic 
investigations of hydrogen bonding have been upon 
liquids, which can show no orientation effects. This 
laboratory, in line with its interest in near infra-red di- 
chroism, has in past years made polarized infra-red studies 
of solid specimens containing oriented hydrogen bonds.'~* 
Such studies are again being undertaken, now with the 
dicarboxylic acids. 

A striking example of dichroism was found in 6-succinic 
acid. Cleavage sections parallel to the crystallographic 
c axis and approximately 0.2 mm thick were used. The 
first spectrograms with polarized light, the H-vector 
vibrating, respectively, parallel and perpendicular to the 
¢ axis, were taken in March 1947. No absorption in the 
1.44 unperturbed O—H region (second harmonic) was 
observed in either run. The run made with the E-vector 
parallel to the c axis was quite unique; strong absorption 
started at roughly 1.57, increased to opacity at 1.87y, 
and remained completely opaque to beyond 2.7y, the 
limit of our recording quartz spectrograph.‘ The per- 
pendicular run showed, in contrast, general transparency 
save for regions of CH2 overtone and combination band 
absorption. These latter were of normal intensity for 
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specimens of the thickness used, the second harmonic CHe 
stretching vibrations, for example, showing about 40 per- 
cent absorption. 

A similar investigation in the fundamental vibration 
region for a high temperature modification of succinic 
acid has recently been reported by D. A. Crooks.*® He found 
a great increase in the absorption of the broad band ex- 
tending from 3250 to 2300 cm ! as the #-vector was rotated 
from the perpendicular to the parallel direction. His con- 
clusion, like ours, is that intermolecular hydrogen bonds 
are responsible for this broad unidirectional absorption. 

The intensity of Crooks’ broad hydrogen bond absorp- 
tion in the fundamental region is, if we interpret correctly, 
less than } that of the fundamental CH: valence stretching 
vibrations at 2900 cm™. A contrast of these findings with 
our broad region of complete opacity indicates an anoma- 
lous ratio of overtone (and combination) to fundamental 
mode intensities for the hydrogen bond vibrations. 

Data from x-ray diffraction® indicate that both the 
chain axis of the succinic acid molecule and the direction 
of the intermolecular O—H---O bonds are very nearly 
parallel to the crystallographic c axis. As was remarked by 
Crooks, the polarized infra-red story bears this out. 

X-ray diffraction does not locate definitely the light 
hydrogen nuclei, their positions on the line connecting the 
oxygen nuclei of neighboring molecules being the reason- 
able inference. This inference is strongly supported by the 
unidirectional nature of the hydrogen bond absorption. 
It is of interest that no such clear-cut confirmation was 
found earlier at this laboratory for the looser hydrogen 
bonds in pentaerythritol.” 

Adipic acid crystals were also examined with results 
similar to but not so marked as those of succinic acid. 
Since the ¢ axis of adipic acid crystals, to which the hy- 
drogen bonds are again parallel, is inclined to the surface’ 
rather than in it as for succinic acid, conditions are not 
so favorable. 

We have been attempting, up to now with small success, 
to grow good crystals of anhydrous a- and #-oxalic acid. 
These, particularly the 8-form, should also show the intense 
hydrogen bonding dichroism. Since there would be no CHz 
absorption in the perpendicular vibration direction, a suit- 
able crystal might be useful as a near infra-red polarizer. 

We attempted to make spectrograms of oxalic acid 
dihydrate specimens but were frustrated by a curious 
phenomenon. Clear sections became very quickly opaque 
in the volume through which the converging beam had 
passed. Presumably water molecules or perhaps ionized 
portions of these were knocked out of the crystal by the 
beam in a sort of photo-molecular effect. 

A detailed account will be published in the future when 
observations and the working out of selection rules from 
the theory of group representations will have been com- 
pleted. A discussion of CHe as well as O—H vibrations 
will be included. 

1J. W. Ellis and J. Bath, J. Chem. Phys. 6, 221 (1938). 

2J. W. Ellis and J. Bath, J. Chem. Phys. 7, 862 (1938). 

3 J. W. Ellis and J. Bath, J. Am. Chem. Soc. 62, 2859 (1940). 

4J. W. Ellis, Rev. Sci. Inst. 4, 123 (1933). 

5D. A. Crooks, Nature 160, No. 4053, 17 (1947). 


6 Llewellyn, Cox, and Goodwin, J. Chem. Soc., p. 883 (1937). 
7W. A. Caspari, J. Chem. Soc., p. 3235 (1928). 
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A Spectroscopic Evidence for Activation of 
Fluorescence by High Valency 
Manganese Ions 


G. SziGeti, E. NaGy, AND E. Makal 
United Incandescent Lamp and Electrical Company, Ltd., 
Vipest 4, Hungary 


October 15, 1947 


E have measured the emission spectra of manganese 

activated zinc beryllium silicate and zinc silicate 
phosphors. The spectrum represented by the number of 
emitted photons as a function of the frequency was found 
to be a sum of several gaussian distribution curves, the 
maxima of each gaussian being at the same frequency for 
both materials (hy=2.35, 2.27, 2.17, and 2.05 electron 
volts), the relative intensities of the individual components 
being, however, different. The main band of zinc silicate 
was found at Ay=2.35 and that of zinc beryllium silicate 
at hy=2.05 electron volts. 

Similarly analyzing the absorption spectra of different 
manganese compounds in solution, we have found the re- 
markable fact that the extinction coefficient of K2MnQ, in 
a strong alkaline solution gives rise to the same curves as 
the emission spectrum of zinc beryllium silicates (see Fig. 
1). 

The analysis of the absorption spectrum of the aqueous 
solution of KMnQO, and Zn/Mn0O,/2 shows the same bands 
as zinc silicate, but in this case the intensity ratio of the 
individual bands is different. The position of the bands 
was found independent of the cation, being determined 
only by the anion. 

These data led us to the conclusion that the emission 
centers of the fluorescence are associated with the (MnO,)~~ 
and (MnQ,)~ ions imbedded in the silicate lattice. This 
means that in the divalent manganese compound added 
for activation there is the necessary quantity of the per- 
manganate resp. manganate ions present, or it will be 
formed during the applied treatment. The number of the 
negative ions necessary for fluorescence can be estimated 
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of zinc beryllium silicate; ordinate scale for both curves gives the same 
maximum, 


7 
6) 
1e ‘ 
se 
of 
is 
id 
X- 
re 
1e 
to 
ar 
e- 
yn 
ne 
Ss, 
). 
ls 
} 
ts 
rd 
ic 
is 
li- 
es 
-3 
ic 
ic 
ne 
or 
he 
as 
or . 
yn 
he 
id 


882 LETTERS TO 


from the saturation value of the emitted fluorescence and 
by assuming the time of the elementary processes. 

For both of these there can be found in the literature 
only very feeble and uncertain data, but they can still be 
relied upon for an approximate calculation that the neces- 
sary number of centers should lie between 10!4+-10" centers 
per gram of fluorescent material. 

This means that if the 10~‘th part of the present man- 
ganese atoms has a valency of six resp. seven, a sufficiently 
great number of manganate resp. permanganate ions is 
present to give rise to the measured fluorescence. It is to be 
noted that the estimated number of these colored anions 
per unit surface does not cause any perceptible discolora- 
tion even in aqueous solutions. Therefore these quantities 
of the activator anions cannot be detected in the powders 
by absorption measurements. 


Raman Effect of Dichlorohexafluorocyclobutane 


WALTER F. EDGELL AND FRANCIS E. KITE 


Department of Chemistry and Chemical Engineering, 
State University of Iowa, Iowa City, Iowa 


November 13, 1947 


ATELY the question of the planarity of the ring in 
simple cyclic hydrocarbons has been of considerable 
interest. Recent work! indicates that the cyclopentane ring 
is puckered. On the other hand, Raman and infra-red 
spectroscopy? indicates that the cyclobutane ring is planar 
while the same seems to be true of the ring in perfluoro- 
_cyclobutane.’ This study of dichlorohexafluorocyclobutane 
was undertaken in connection with this laboratory’s interest 
in this problem. 

Dichlorohexafluorocyclobutane is a colorless liquid boil- 
ing at 60°C. The material used was furnished by the E. I. 
DuPont de Nemours and Company. It was further purified 
by vacuum distillation, dried over CaCl for twenty-four 
hours, and made optically clear by shaking with activated 
charcoal and filtration through a very fine sintered glass 
filter. The Raman tube, having a volume of 15 ml, was 
filled by vacuum distillation. 

Exposures were made on Eastman spectrographic plates 
type 103-0 and type 103a—0, using an Adam Hilger E-2 
spectrograph. The Raman tube was suspended in a vertical 
position with the scattered light being directed into the 


TABLE I, Raman lines of dichlorohexafluorocyclobutane. 


Estimated Estimated 
Ap relative Av relative 
intensity intensity 
2 

162 Ss 863 W 
216 VW 910 W 
270(b) 966 Ss 
334(b) Ss 1074 W 
396(d) VS 1111 WwW 
465 S 1190 W 
490 S 1243 M 
561 Ww 1287 W 
640(b) VS 1392 M 
810 VW 


S =strong; M =medium; W =weak; V =very; b =broad; d =diffuse. 
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spectrograph by a right-angle prism. Eight upright mercury 
arcs arranged in a circle about the Raman tube provided 
the exciting radiation. Each of these mercury arcs were 
mounted in front of an individual elliptical reflector which 
was focused upon the Raman tube. The 4358A line was 
used as the exciting intensity. Inasmuch as fluorocarbons 
are notoriously poor light scatterers, exposure times varied 
from 48 to 172 hours by use of slit widths from 25 to 60 
microns. A continuously circulating filter solution of 0.2 
percent paranitrotoluene and 0.03 gram of DuPont Rhoda- 
mine 5 GDN Extra per liter, dissolved in alcohol, was used* 
to reduce general scattering in the region of 4358A and to 
eliminate the Raman shifts excited:-by the 4047A mercury 
line. The average temperature of the sample was 38°C 
during the exposures. 

We have observed nineteen Raman shifts in dichloro- 
hexafluorocyclobutane. The results are shown in Table I. 

We wish to thank Dr. J. D. Park of the Jackson Labora- 
tories of the E. I. DuPont de Nemours and Company for 
this substance. Further consideration of these results will 
be published later. 


1J. E. Kilpatrick, K. S. Pitzer, and R. Spitzer, J. Am. Chem. Soc 
69, 2483 (1947); J. G. Aston, S. C. Schumann, H. L. Fink, and P. M. 
Doty, J. Am. Chem. Soc. 63, 2029 (1941); J. G. Aston, H. L. Fink, 
and S. C. Schumann, J. Am. Chem. Soc. 65, 341 (1943); E. J. Rosen- 
baum and H. F. Jacobson, J. Am. Chem. Soc. 63, 2841 (1941). 

2T. P. Wilson, J. Chem. Phys. 11, 369 (1943). 

3W. F. Edgell, J. Am. Chem. Soc. 69, 660 (1947). 

4J. T. Edsall and H. Scheinberg, J. Chem. Phys. 8, 520 (1940). 


Phases of Fourier Coefficients Directly 
from Crystal Diffraction Data 


D. HARKER AND J. S. KASPER 


Research Laboratory, General Electric Company, 
Schenectady, New York 


July 14, 1947 


HE intensity of an x-ray beam diffracted from a 
crystal can be used to find the absolute value of the 
“crystal-structure factor” | F,x:|, but not its phase. How- 
ever, it appears to be possible to place limits on the phase 
or, if Fyx: is real, to determine the sign, by combining 
the intensity data from several diffracted beams. The 
general methods for doing this are outlined in the following 
paragraphs; a detailed article on the subject will appear 
later. 

Although not necessary, it is helpful in strengthening 
the method to assume that the “‘atomic-structure factors,” 
fj, of the atoms in a crystal are related to one another by 
the relation f;=Z;f, where Z; is the atomic number of the 
j* atom and f is a function of (sin @)/d called “the unitary 
atomic-structure factor.” If there are N atoms in the unit 
cell of a crystal and 


is the total number of electrons in a unit cell, then 


f;=Zf. 


j=l 


‘ 
— 
_ 
i, 
N 
Z=2>7Z; 
= 
{ 
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Since 


N 
Frei = 2 fi exp [ 
the introduction of the f assumption allows one to write 

N 

= nj exp [ 

where Fet=(Fare)/(Zf) and n;=Z;/Z. Par is called “the 
unitary crystal-structure factor’ or “the unitary Fourier 
coefficient.” This way of writing assumes that all the atoms 
in the crystal have the same “‘shape.”’ 
The relation 


N N N 
j=l j=l j=l 


(“‘Cauchy’s inequality’’) can be applied to the expression 
for Fix. Taking a; = (n;)t and b; = (n;)} exp [ —2mi(hx;+ky; 
+/z;)], one obtains 


N N 
| nyt exp [—2ri(2hx } 
j= 


or | since 


This result explains the use of the word “unitary” in the 
name of 

If the crystal has symmetry, more interesting relations 
are obtained from Cauchy’s inequality. If, for instance, a 
center of inversion is present, one can write 


2 nj cos 


i>! 


and, noting that Fy. is real and applying Cauchy’s in- 
equality, 


N/2 N/2 
i=l i=l 


By the use of 2 cos?a=1+cos 2a and elementary rear- 
rangements, it follows that This im- 
portant relation requires, for instance, that Fej2%2. be 
positive or zero if or, if | then 
must be positive if F%,,..2=4, ete. Conversely, if this in- 
equality is contradicted, the crystal cannot have a center of 
symmetry. (This does not invalidate Friedel’s Law, which 
merely states that | Fixi| = | Fixi|.) 

The presence of other symmetry elements leads to other 
inequalities which are derived from the appropriate ex- 
pressions for Fixx: by using Cauchy's inequality as in the 
last paragraph. Inequalities for crystals having the several 


TABLE I. The inequalities arising from the various axes of symmetry. 


Symbol 
of Equivalent 
Axis Symbol Coordinates 


Inequality 


142K -0\|cos 2maH-K H+2K-0+4 |\Pu-K H42k.2L|cos 2maH-K H42K. 2L* 


L 
cos H2K-o+} cos HK-O 


(—1)4 


1 Triclinic | S1 
2 — Monoclinic | Prat]? <4 +3 x02 
2 m Monoclinic | Pre]? <3 +32 oxo 

21 Monoclinic | 
3 Hexagonal n42k.0|cos 2man-K H42K-0* 

4 Tetragonal | <4 +4 | cos nok 

(—1)? l 
41 Tetragonal | Pani]? +4 cos Pr 
(-1)! 
3 
6 _ Hexagonal 
m 

61 _ Hexagonal 

L 
62 Hexagonal |FHK. cos H+2K. 0+} cos 2a 

63 Hexagonal H42K.0+ 


6 


Fux.o 


* Note: may be expressed exp is the “phase” of 
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rotation axes, screw axes and rotation-inversion axes are 
listed in Table I. 

It is possible to obtain information on the signs of those 
Fyx1’s not determined by the inequalities involving only 
one | F,,.|%. For instance, for the case of a center of inver- 
sion, treatment of | Pax Fix |? by the methods just de- 
scribed yields the relation 
— Fy | S144 + — — Pee. 
It is probable that further development of this method 
will provide sign or phase determinations for most of the 
Fixw’s of a crystal. That not all of the signs (or phases) 
can be determined by these methods follows from the 
ambiguities in crystal structure determination discussed 
by A. L. Patterson.! 

It is essential, before using any of these inequalities, 
that the experimental F,:’s be expressed in absolute units. 
In the case of x-ray diffraction, the units are electrons per 
unit cell. 


1A. L. Patterson, ‘‘Ambiguities in the x-ray analysis of crystal struc- 
tures,"’ Phys. Rev. 65, 195 (1944). 


Infra-Red Dichroism in Aligned Polythene 
and “Parowax”’ 
LEONARD GLATT AND JOSEPH W. ELLIS 


Physics Department, University of California, Los Angeles, California 
November 6, 1947 


RIENTED specimens of polythene were examined 

for dichroism in the 0.6- to 2.7-4 region. These 
specimens had been very kindly given to us by Dr. Wilfried 
Heller, now at Wayne University, and his former co- 
worker at the University of Chicago, Hans Oppenheimer. 
They were plastic films which had been stretched until 
the alignment of the molecular chain axes was as closely 
parallel to the stretch direction! as possible. The planes of 
the carbon backbones had a random angular distribution 
about this direction. 

The absorption bands which appeared with unpolarized 
light were as follows: (1) a weak band at 8237 cm™; 
(2) a well defined doublet at 5782 and 5671 cm; (3) a 
very strong doublet at 4322 and 4271 cm™; (4) a com- 
paratively weak but distinct band at 4216 cm™; (5), a 
shallow poorly defined region between roughly 4180 and 
4080 

The most interesting dichroism was exhibited by the 
4216-cm™ band. All the bands were present to some extent 
for both E-components, vibrating respectively perpen- 
dicular and parallel to the stretch direction, but were 
considerably stronger in the former. The 4216-cm™ band, 
in contrast, appeared only for the E-vector parallel to the 
stretch direction. Its intensity was considerably stronger 
than that found with unpolarized light. 

Since this band appears distinctly for films as thin as 
50u, it probably involves displacements within the CH: 
groups. It is quite likely an overtone of a fundamental 
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mode in which the change of dipole moment is completely 
parallel to the chain axis. In their mathematical treatment 
of an infinite >CHez chain, Whitcomb, Nielsen, and 
Thomas? derive a fundamental mode in which all the 
hydrogens vibrate against all the carbons in a direction 
parallel to the chain axis. They attribute, we believe 
mistakenly, the 1460-cm™ hydrocarbon band to this mode. 

The known strong bands of normal saturated hydro- 
carbons at approximately 2925, 2853, and 1460 cm™ are 
usually assigned, on the uncoupled oscillator approxima- 
tion, to internal vibrations of the >CHp units. They are 
considered to be, respectively, the va, vs, and 6 type modes 
of a triangular YX» molecule. On this same naive ap- 
proximation, the two pairs of doublets (2) and (3) can be 
readily assigned to 2va, 2vs, va +6, and v.+6 in that order. 
In keeping with these assignments, the high frequency 
component of both doublets was considerably the stronger 
of the two. 

The CHs groups lie on planes perpendicular to the 
chain axis. For uncoupled oscillators it would be expected 
that harmonics and combinations of va, vs, and 6 type CH» 
vibrations would be inactive for the E-vector parallel to 
the chain axis. Although appreciably diminished, both 
doublets definitely appeared for this E-component. Their 
intensity was greater than what can be accounted for by 
imperfect alignment of the polarizer and the fact that we 
used a converging beam. A reasonable explanation is that 
for coupled oscillators there is a component of these vibra- 
tions perpendicular to the CH: planes. It is of interest that 
in both pairs of doublets the lower frequency (presumably 
“symmetrical” ) component showed considerably the lesser 
diminution for E-vibrations parallel to the stretch direction. 

Unpolarized light spectrograms of thin films sliced from 
a block of ‘“‘parowax” were essentially identical with those 
for polythene. The random molecular orientation in the 
sliced films made useless any polarization study. 

Electron diffraction studies of very thin paraffin films, 
deposited from a drop of ether solution onto a water sur- 
face, have shown that the molecular chain axes are oriented 
normal to the surface.? We prepared semi-crystalline films 
of from 40- to 500-u thickness by slow cooling of a melt 
ona hot water surface. To minimize disorientation owing to 
the heat of the beam, these films were placed at the exit 
rather than entrance slit of our instrument. Unpolarized 
light spectrograms gave convincing proof of an excellent 
alignment of the chain axes normal to the surface. The 
band at 4216 cm™ was completely missing, and the in- 
tensity ratios and general appearance of the observed 
bands were unmistakably like those of the spectrograms 
made for polythene when the E-vector was directed per- 
pendicular to the stretch direction. 

A detailed account will be published in the future. We 
intend to investigate the fundamental region with an in- 
strument now under repair. An assignment in agreement 
with group theory selection rules will be attempted for the 
bands involving CH: vibrations. 


1W. Heller, Phys. Rev. 69, 53 (1946). 

2S. E. Whitcomb, H. H. Nielsen, and L. H. Thomas, J. Chem. Phys. 
8, 143 (1940). 

3 Trillat and Hirsch, Comptes rendus 195, 215 (1932). 
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On the Frequency Factor in the Viscosity- 
Temperature Relationship of Liquids 


M. S. TELANG 


Laxminarayan Institute of Technology, Nagpur University, 
Nagpur, India 


November 5, 1947 


N recent publications'? I clarified the term ‘activation 
energy” for viscous flow of liquids from the expression 


n=Z exp [y(M/D—d)}- N1/(cos30 °)#/RT]. (1) 


Since it is possible to calculate Eyis, directly according to 
(1), it is necessary to examine whether it is also possible 
to calculate the “frequency factor,” Z, directly. Eyring 
et al.? have derived the equation 


= (F/Fy)-exp Leo/kT], (2) 


the symbols being the same as the original ones of Eyring 
et al. The product \2AsA;1 has been assumed to be approxi- 
mately the volume inhabited by a single molecule in the 
liquid state and, hence, it has been put equal to V/N, 
where V is the molar volume. Eyring ef al.* have further 
transformed Eq. (2) into 


n=(hN/V)-expAEvap,/2.45RT (3) 


Now, I find that the non-exponential term representing 
the “frequency factor’’ has to be also altered. The modified 
equation obtained by me is 


n=(hN/V—b)-exp [y(M/D —d)}- N3/(cos 30°)!/RT], (4) 


introducing a constant b, characteristic of each substance. 
The symbols in (4) are:—y=viscosity, y =surface tension, 
M~=theoretical molecular weight, D=density of liquid, 
d=density of vapor, N=Avogadro number, R=gas con- 
stant (in ergs), T=absolute temperature, Planck’s 
constant, and V=molar volume. The only unknown factor 
b in (4) can be easily evaluated by making appropriate 
substitutions for the experimental values of 7, V, 7, D, d, 
and 7, at one temperature. Calculating b from (4) in this 
manner, Table I shows that b/ V; is nearly the same for all 
the substances listed for illustration. The ratio of b to V., 
the molar critical volume, is very nearly }. It follows, 
therefore, that the factor 6 introduced in (4) is almost 
identical with van der Waals’ constant b. This observation 
applies to other non-associated substances also and is 


valid at temperatures even far above the melting points. 


The necessary data have been taken from the Inter- 
national Critical Tables and from Lange’s Handbook of 
Chemistry.® 

The correction factor 6 introduced indicates that A2sA1 
has to be interpreted as the free volume or a hole in the 
liquid and not as the total volume inhabited by a molecule. 
It is quite likely that hN/(V—b) may be a mathematically 
simplified form of 


V—b —b 
The complete equation would then be 
hN b 
1+ exp [y(7/D —d)}- N4/(cos 30°)/RT]. 


(S) 
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TABLE I, 

Substance b/Ve 
0.313 
Benzen 0.318 
Coben 0.338 
Chloroform 0.338 
Cyclohexane 0.346 
Ethylene bromide 0.314* 
Ethyl ether 0.344 
n-Hexane 0.336 

0.331 


* V. is calculated as thrice van der Waals’ constant. 


A somewhat similar idea has been incorporated by Eyring 
et al.® in giving an alternative treatment of viscosity by 
the theory of absolute reaction rates. It is rather premature 
at this stage to arrive at any definite conclusion regarding 
the exact significance of the correction factor. However, a 
mathematical analysis of the factor may be briefly given 
here. The number of holes , per mole of liquid is equal to 
V—b/vn, where v is the increase in volume required to 
form a single hole.’ If b/N is the volume of a single mole- 
cule, then 7,=6/N; since Eyring’s interpretation of the 
“activation energy” as AEyap./n=Evis. is no longer ap- 
plicable, there seems to be no justification for assuming 
that v, is not equal to b/N. Hence, 


V—b/b-N. 


b 
( nh 


1M.S. Telang, J. Chem. Phys. 15, 525 Had 

2M. S. Telang, J. Chem. Phys. 15, 844 (1947) 

2S. Glasstone, K. J. Laidler, and H. Eyring, The Theory of Rate 
Processes, (McGraw-Hill _ Company, Inc., New York, 1941), p. 483. 

4 See reference 3, p. 493 

5N. A. Lange, Handbook of Chemistry, (Handbook Publishers, Inc., 
Ohio, 1944). 

6 See reference 3, p. 488. 

7 See reference 3, p. 510. 
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Prism Spectrometry from 24 to 37 Microns 


EARLE K. PLYLER 
National Bureau of Standards, Washington, D. C. 
November 14, 1947 


Y the use of a prism of KRS-5 (thallium bromide- 

iodide) with a refracting angle of 26°, it is possible to 
make measurements in the infra-red to 37 microns. With 
a thinner prism it may be possible to extend the observa- 
tions to 40 microns. The prism was incorporated in a 
Perkin-Elmer spectrometer. The stray radiation is almost 
completely removed by the use of reflection filters of LiF 
and CaF». Greater deflections will be obtained by using a 
filter of NaF beyond 36 microns. The instrument was 
calibrated by the use of known absorption bands in this 
region. The dispersion of KRS-5 is very high from 0.546 
to 1.00 micron, but from 2 to 15 microns it becomes ex- 
tremely small. The dispersion is very good from 24 to 37 
microns. 
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POLYSTYRENE 


POLYETHYLENE 
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WAVELENGTH IN MICRONS 


Fic. 1. Transmission of polystyrene 25 microns thick, polyethylene 
240 microns, KRS-5 4 mm, KBr 6.4 mm, and AgCl 2.1 mm. 


In Fig. 1 is shown the transmission of several materials. 
The sample of KRS-5 is 4 mm thick and transmits about 
70 percent at 24.5 microns. This is about the amount trans- 
mitted at 20 microns, and most of the loss is due to reflec- 
tion. As the wave-length is increased, the transmission 
decreases, and is 45 percent at 37 microns. Some samples 
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of KRS-5 have less transmission than the sample shown 
in Fig. 1. One sample examined is opaque beyond 32 
microns. This difference between samples obtained from 
different sources may arise from impurities in the crystals. 
Polystyrene shows the highest transmission in this region 
of any substance measured. This is partly caused by the 
fact that the layer is only 25 microns thick. When a 1-mm 
thick layer is used only 25 percent of the incident radiation 
is transmitted at 33 microns. However, for coatings and 
windows thin films of polystyrene are useful. A layer of 
polyethylene 240 microns in thickness also shows high 
transmission in this region. Thinner layers show somewhat 
higher transmission. The irregular form of the transmission 
curve is probably caused by absorption bands in the ma- 
terial. Films of C2F4, CoH2F2, and C2H2Cle were also meas- 
ured but they have low transmission in this region. The 
KBr sample is 6.4 mm thick and transmits about 10 percent 
at 31.5 microns. Thus, in thicknesses of one mm or less 
KBr will be useful as a window to 30 microns. A sheet of 
AgCl, 2.1 mm thick, transmits 10 percent at 27.5 microns. 
However, the sample is slightly discolored and may not 
have as high transmission as would be found for a clear 
specimen. Now that the instrument has been calibrated, 
investigations are being carried out on a number of liquids 
and vapors. 
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